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FOREWORD 


The computer program presented in this report was developed for the 
Vibration and Acoustics Branch of the George C. Marshall Space Flight 
Center, National Aeronautics and Space Administration. This work has been 
accomplished by the Vibration and Acoustics Group, Structural Engineering 
Branch, Chrysler Corporation Space Division, Huntsville Operations. 



SUMMARY 


This computer program is to calculate the random vibrational response 
of a rectangular cylindrical shell panel cross-reinforced with ribs and 
stringers subjected to the excitation of fluctuating pressure environments* 

Three boundary conditions are considered* four edges simply supported; four 
edges clamped; and two opposite edges simply supported while the other two 
clamped. The special cases of a complete cylinder and a flat panel are 
included. This program is written so that either all three boundary conditions 
or any one boundary condition can be selected m any run. The normal mode 
approach is used m the formulations. Responses calculated are the acceleration, 
displacement, and stress spectral densities. Mean-square and root-mean-square 
values are also calculated. Output data are tabulated and plotted. This 
Manual is written according to the documentation requirements specified by 
the Computation Laboratory of MSFC. It contains three sections £es cnbxng the 
problem, the programming, and the deck setup m detail. 
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INTRODUCTION 


Recent advancements m the mathematical formulations of solutions for 
structural responses due to random loading has made it feasible to develop 
computer programs for use by the structural engineer* The computer program 
reported in this manual calculates the random vibrational responses of rec- 
tangular cylindrical shell panels cross-reinforced with ribs and stringers* 

It is the result of a research project* 

The program is written m FORTRAN IV language for the IBM 7094 computer* 
The results of this project are reported m two volumes. 

1. Volume I - Theoretical Analyses 

2* Volume II - User’s Manual 

The first volume contains the equations, technical discussion of the 
program, and analyses of the results* The second volume is oriented to utili- 
zation of the programs. 
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SECTION I. PROBLEM 


A. Abstract 


This is a program to compute the vibrational responses of a rectangular 
cylindrical shell panel cross-reinforced with stifeners, subjected to the 
excitation of fluctuating pressure environments. The boundary conditions 
considered are: four edges simply supported; four edges clamped, and two 

opposite edges simply supported while the other two clamped. The special 
cases of a complete cylinder and a flat panel are included The normal mode 
is used m the formulations. The responses are calculated as acceleration, 
displacement, and stress spectral densities, and the overall mean-square 
and root-mean-square values. The output data are tabulated and plotted. 

The program is written so that either all three boundary conditions or any 
one boundary condition can be selected m any run. 
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B. TECHNICAL DESCRIPTION 


This is a program written to calculate the random vibrational responses 
of rectangular cylindrical shell panels cross-reinf orced with stiffeners. 

Three boundary conditions are considered: four edges simply supported; four 

edges clamped; and two edges simply supported while the other two edges 
clamped. Special cases of flat panel and complete cylinder are included. The 
one- third-octave spectrum of the excitation pressure is read in as input in 
any discrete frequency. The program will apply when either the whole panel 
or a portion of the panel is exposed to the excitation The normal mode ap- 
proach is used m the formulation. Analytical expressions for the joint 
acceptance are derived for all mode combinations. Contributions of the mam 
terms as well as the cross terms are accounted fof to obtain the responses. 

The responses at any point are calculated as displacement, acceleration and 
stress spectral densities. Mean-square and root-mean-square values of the 
responses are calculated by numerical integration. The response spectral 
densities are tabulated and plotted. The frequency range is from 5 (or lower) 
up to 5000 Hertz. One-nth-octave band is used for the frequency increment. 

The number of data points can be increased up to 500 for each spectral density 
plot. For each data point, 625 terms are summed to give the spectral density. 

In addition to the excitation spectrum, the input data includes the 
geometric dimensions, material properties of the panel, and some control 
constants . 

The program is written so that either all fhree boundary conditions or 
any one boundary condition can be calculated m any run. 

Computed responses have been compared with experimental results, and 
the agreement is very good. 
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C. EQUATIONS 


a* Subroutine RSR - Calculation of Response of Simply- Sup ported Rectangular 

Shell Panels Cross-Reinforced With Stiffeners 

1 . Natural Frequency 

The undamped natural frequencies are given by 
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Refer to Figure 1 for geometric dimensions. 

a = radius of shell panel (input) 

E = Young's modulus of panel (input) 

E' = Young's modulus of stiffeners (input) 

h = Thickness of 'panel skin (input) 

h’ = Smeared-out thickness of stiffeners (input) 

v = Poisson’s ratio of panel (input) 

= Spacing of y-direction stiffeners (input) 

= Spacing of x-direction stiffeners (input) 

1^ = Moment of inertia of one x-direction stiffener with respect 

to neutral axis of cross-section of panel (input) 
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= Moment of inertia of one y-direction stiffener with respect 
to neutral axis of cross-section of panel (input) 

p f “ Mass density of stiffeners (input) 

p = Mass density of panel (input) 

\ 

2. Responses 

The displacement response spectral density at point r(x,y) is given by 
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Circumferential width of panel subjected to excitation (ipput) 
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to = 
f = 
j»k,m,n = 

“jk = 

The magnitudes 


Frequency xn rad/sec, - 2?rf 
Frequency in Hertz (independent variable) 

Mode indices 

Undamped natural frequncies given by Equation (1) 
of the frequency response functions are: 
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The joint acceptances squared for different mode combinations are given by 
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where 


A, = Decay constant m x-directxon (input) 

A = Decay constant in y-direction (input) 

Z * 4 
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The excitation spectral density a,n 


(psi) 

rad/sec 


$ pp (fcj), is given by 


Sp p (f) - 170.576 


W*° ' h (10) 


The excitation spectral density m decibels per Hertz, S (f ) , is given by 

r r 


Spp(f) = S 3rd (f) - 10 log 10 (0.23157f) 


S 3rd (f) 


One-third octave excitation pressure level m decibels 
(input) 


The acceleration response spectral density at r(x,y) m 
v sec ' 


rad/sec 


is given by 
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p- 

The acceleration response spectral density m is given by 


S w (r,f) = 4.215093 x 10 5 * w (£,<o) 


(4) 


The mean-square acceleration m rl g^ rl is given by 
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The root-mean-square acceleration in l, g is given by 


-> / o 
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The stress spectral density m (psi) is given by 

rad per sec 
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largest height of stiffeners at r, see Figure 1 
(input) 
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I [(f) 2 + ’(f) 2 ] 


mirx __ rmy 
sj-n — 


=1,3.. 
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The stress spectral density m (psi) / Hertz is given by 


S art< r » f > 


aTT^-rtCrjOj) 


The mean-square stress is given by 


9 -> 

tf 2 (r) = 


/ S aa^ £)df 


The root-mean-square stress is given by 


cf (r) ~ [ cr 2 (r) ] ^ 


(inch) 

The displacement spectral density m is 

Hertz 


S (r,f) = 2 tt$ (r,w) 
w ww s 9 ' 


The mean-square displacement is 


w 2 (r) = / S (r,£)df 


The root-mean- square displacement is 


w(r) = w^(r) 
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b . SUBROUTINE RFR 


DYNAMIC RESPONSE OF FOUR SIDES CLAMPED RECTANGULAR 
SHELL PANELS CROSS-REINFORCED WITH STIFFENERS 




11 



2 " 1/2 

to . = ir M 

jk 




4 


+ 


2H 


(j + |) (k + I)[(j + |) - 0-6366] [(k +k) - 0. $366] 


A 2 


2 4 
a ir 


1 + 


j,k = 2,3,. 


Eh 

k + A 

2 

3 + 1 


T 2 


(t) 


1/2 


X 


Eh 3 . E,I 1 
+ 


12(l-v)' 


Eh 


12 (1— v j 


E'l„ 


H 


M 


Eh" 


12(l-v ) 
ph + p 'h' 


The natural frequencies m Hertz are given by 



2ir 


E = Young’s modulus of panel (input) 

E’ = Young's modulus of stiffeners (input) 

h - Thickness of panel skin (input) 

h' = Smeared-out thickness of stiffeners (input) 


(ID) 


(IE) 


( 2 ) 
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V 


a 

a 

b 

I 

I 


1 

1 

1 

2 


P 

P 1 


Poisson’s ratio of panel (input) 

Radius of shell (input) 

Spacing of y-direction stiffeners (input) 

Spacing of x^direction stiffeners (input) 

Moment of inertia of one x-direction stiffener with 
respect to neutral axis of cross-section of panel (input) 

Moment of inertia of one y-direction stiffener with respect to 
neutral axis of cross-section of panel (input) 

Mass density of panel skin (input) 

Mass density of stiffeners (input) 


2 „ Responses 


The displacement response spectral density at point r(x,y) is given by 
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(r«w) 


where 


$ (r 5 w) 

ww v 5 ' 


- S ,2 *pp(“> l F lt <J)F mn (?)|H 


H_ 


jk v ' uin v ' * jk 1 1 mn 1 jkmn 


J ,k,m,n 

=1,2,3. . . 


= displacement spectral density in 


inch 


rad/sec 


(3) 


S' = Area of panel subjected to excitation (mch^) 

= b'S.' 

b' = Width of panel subjected to excitation (input) 
l' = Length of panel subjected to excitation (input) 

Fjk(r) = = normal mode 

F mn (r) = X m Y n = normal mode 
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Undamped natural frequencies given by Equation (1) 
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The magnitudes of the frequency response functions and are given 
by equations (2F) and (2G) m Section a. 

The joint acceptance squared for different mode combinations are given by 
equations (2H) through (2L) m Section a. 


The displacement spectral density m (inch) 

Hertz 

S w (r, f ) = 27 r$ ww (r,w) 

The mean square displacement m (inch) 2 is 

w 2 (r) = / S ww (r ’ f > df 

The root-mean square displacement in inches is 
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The excitation spectral density m decibels per Hertz, S n (f ) , is given by 

FP 

s pp (f) - s 3 rd (f) ■ 10 lo S 10 (0-23157f) 

where 


S 3 rd(f) - One- third octave excitation pressure level m 

decibels (input) 

< 2 

The excitation spectral density in (psi) /Hertz is given by 
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The acceleration response spectral density at r(x,y) in 

\ sec 2 / rad/sec 


is given by 


\<fW 




The acceleration response spectral density m B 1S given by 

Hertz 


s ww( r » f ) 


4.215093 x 10 5 $^( 1 , 11 )) 


o 

The mean-square acceleration m !, g " is given by 
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The root-mean-square acceleration in "g" is given by 
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The stress spectral density m Cpsi) is given by 
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where ^(r,w) is given by Equation (3) and y Z (r) is given by Equation (7A) 
m Section a. 


The stress spectral density in 
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The mean-square stress m (psi) Z is given by 
c 2 (r) = / S o0 (r,f)df 
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The root-mean— square stress m psi is given by 


o(r) 
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( 19 ) 



C. SUBROUTINE RSF 


DYNAMIC RESPONSE OF TWO OPPOSITE EDGES SIMPLY-SUPPORTED AND OTHER 
TOO CLAMPED RECTANGULAR SHELL PANELS UNDER RANDOM PRESSURE FIELD 

Simply-Supported Edges: x = 0, x = £ 

Fixed Sides : y = 0 , y = b 

1. Natural Frequencies 

The natural frequencies m radians per second are given by 
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The natural frequencies m Hertz are given by 


(IE) 



2ir 


( 2 ) 


E = Young's modulus of panel (input) 

E' Young's modulus of stiffeners (input) 

h = Thickness of panel skin (input) 

h' = Smeared-out thickness of stiffeners (input) 

v = Poisson's ratio of panel (input) 

a = Radius of shell (Input) 

a^ = Spacing of y-direction stiffeners (input) 

= Spacing of x-direction stiffeners (input) 

I^ = Moment of inertia of one x-direction stiffener with respect 
to neutral axis of cross section of panel (input) 

I 2 = Moment of inertia of one y-direction stiffener with respect 
to neutral axis of cross section of panel (input) 

p' = Mass density of stiffeners (input) 

p = Mass density of panel skin (input) 


2. Responses 

The formulation for responses for Subroutine RSF are the same as Subroutine 
RFR given m Section b, except the normal modes are by the following expressions; 
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D. DEFINITION OF TERMS 

MNEMONICS 

FORMULAS 

DESCRIPTION 

SPPF(I) 

^pp ) 

Excitation spectral density m decibels per 
Hertz. 

FIPW(I) 

$ pp <“> 

2 

Excitation spectral density in (psi) rad per 
sec. 

FNW(J,K) 

w Jk 

Nautral frequencies of panel 

OMEGA 

0) 

Independent frequency variable for spectrum 
m rad/sec 

POWJ2 

J 2 

jkmn 

Joint acceptance squared 

PIWW 

$ (r.w) 

ww 5 ' 

n 

Displacement spectral density in mch z /rad 
per sec 

OMEG(I) 

f 

_v 

Independent frequency variable for spectrum 
m Hertz 

SWW(I) 

“T 

S (r,f ) 
ww 

2 

Displacement spectral density m m /Hertz 

PIWG 

ww v * 1 

2 4 

Acceleration spectral density m m / sec / 
rad per sec 

PIWGl(I) 

-¥ 

O 

Acceleration spectral density in g^/Hertz 

PSSW 

^cjaCr.to) 

2 

Stress spectral density in psi /rad per sec 

SSSF(I) 

Sera (r,f ) 

2 

Stress spectral density m psi /Hertz 

SPPP (I) 

s pp«) 

2 

Excitation spectral density m psi /Hertz 

QX 


2 

Quantity for the calculation of y 

QY 

Q y 

2 

Quantity for the calculation of y 

QW 

Qw 

2 

Quantity for the calculation of y 

GAM2 

Y 2 (r) 

Constant to change displacement spectral 
density into stress 

PAS 

2 

P 

a 

Overall mean square excitation pressure 

2 "*■ 
W (r) 

ATS1 

Mean square displacement 
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D. DEFINITION OF TERMS (Continued) 


MNEMONICS 

FORMULAS 

DESCRIPTION 

a 2 (r) 

ATS 2 

Mean square stress 

g 2 (?) 

ATS 3 

Mean square acceleration 

W(r) 

ATI 

Root-mean square displacement 

- 3 - 

a(r) 

AT2 

Root-mean square stress 

G(r) 

AT3 

Root-mean square acceleration 
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E. SPECIAL OPTIONS 


This program can be easily modified to caluclate the acceleration 
spectral density in decible scale referenced gravity acceleration, and 
compute the vibro-acoustic transfer function Modification to calculate 
the average responses over the whole structure and to investigate the contri- 
bution of the cross terms to the response is also not difficult, 

F. NUMERICAL METHODS OF SOLUTION 

The calculation of the mean-square responses is by numerical integration 
of the area under the spectral density curve. However, the integration pro- 
cedure is written inside the program and no integration subroutine is required, 

G. TECHNICAL REFERENCES 

1 Lee, T. N. , "Computer Program for Prediction of Structural Vibrations 
to Fluctuating Pressure Environments, 1 ' Contract No. NAS 8-21403, 

Monthly Progress Reports 1 Through 11. August 1968 to June 1969. 

H. RELATED PROJECTS 

This computer program has been used m the project, "Comparative Analysis 
of Acoustic Testing Techniques, 11 MSFC Contract NAS8-21425, which Chrysler 
Huntsville Operations is conducting from July 1, 1968, to July 31, 1969. Com- 
parison of the computed responses with the experimental data shows good 
agreement . 
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SECTION II . PROGRAMMING 


A. Library Subroutines 

No non-system subroutine is required for this program, 

B Program Subroutines 

This program contains a short mam driver program and five subroutines 
which are called xn the following order: 

1 Subroutine PENT is used to print the output data. 

2. Subroutine GRIDIV is a general grid subroutine for the SC 4020 
and is called by the plot routine in subroutines RSR, RFR and RSF 

3. Subroutine RSR is to calculate the frequencies and the responses 
of simply supported rectangular panels. 

4. Subroutine RFR is to calculate the frequencies and the responses 
of the four -"edges -clamped rectangular shell panels. 

5. Subroutine RSF is to calculate the frequencies and the responses 
of panels with two edges simply supported, the other two clamped. 

C. Special Input Tapes * 

None 

D Special Output Tapes 

A-8 is the plot tape used by the SC 4020 plotter. 

E. Plots Generated 

The following eleven plots are generated by this program. In all the 
plots, the frequency is the independent variable while the spectral density 
is the dependent variable 

a Spectral density of the excitation pressure in decibles versus 
frequency xn Hertz. 

2 

b. Spectral density of the excitation pressure m (psi)^ /Hertz 
versus the frequency m Hertz. 

The following three plots are generated for each of the three boundary 
conditions 

c. Acceleration spectral density m g^/Hertz versus frequency m Hertz. 

d. Displacement spectral density m mch^/Hertz versus frequency m Hertz. 

e. Stress spectral density in (psi) /Hertz versus frequency in Hertz. 


24 



F. BLOCK DIAGRAMS 
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c. SUBROUTINE RSR 














d . SUBROUTINE RFR 















e. SUBROUTINE RSF 















G. FLOW CHARTS 
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b. SUBROUTINE GRIDV 
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c. SUBROUTINE RSR 


/sUBROliTINL 

V nsa 


0 


INITIALIZE 
BCDX.BCDY, 
BODY 1, BCDY2, 
AND BCDYM 
WHITE HEADING 
R* WIND 3 


CM I CAMUAV 
CM LSMXYY 
C\I LSCOUTY 
SI T CONST AN IX 
ICT,JK AND PI 


9 


CYLCULATE 
SPPF (I), 
FIPW (I) 



COMP 

cons: 

HC, D^ 
AND I 

UTE 

rANTS 

k^DX 

m 


20 


COM 

FHZ 

FNW 

PLTF 

AND 

(h3) 

49 

: 

COMPUTE 
S,QX QY 
AND QW 


COMPUTE HI, 
GAM 2, OMEGA 
1CT = 0, IC = 1 


BEG 

LOO 

IN 21 
P 



BEG 

LOO 

IN 40 
P 



KA = JA 

COMPUTE 

PIPP 

4-> ■ 


r 

[ PIWW - 0 0 



BEGIN 22 
LOOP 




COMPUTE 
FM1, FM2, 
RMA, HJK, 
HMN, AC, 
BC, CC, 
POM EG, 
POWJ2, 
PIWW 



CONTINUE 

a 



CALCULATE 
OMEG (I) 
SWW (I) 

PIWG 1 (I) 
SSSF (I) 

SPPP (I) 



IC = 2 

FT - OMEG (I) 
G1 = SWW (I) 

G2 = SSSF (I) 

G3 - PIWGl (I) 
ATS1 = 00 
ATS2 = 00 

ATS3 - 0 0 




(Am la 1 I 

A 1 SI, AISJ, 

A IS 

i, Ail, 

AT 2 

\ I i 

PI 

OMf (,(1) 

G1 

SWW [If 

02 

SSSI (i) 

OI 

piwc.ki) 


CHICK I OR 
MIN-M \\ 

Y All I S OI 
SWW (I, 

SSSI (1, 
PIWOlUi 






com ini i 

KM I 



98 ^ 



END FILE 3 
REWIND 3 
CALL PRNT 
SET PLOT 
SCALES 


301 


CALL SMXYV 
PLOT SPPF, 
SWW, SSSF, 
SPPP, PILG1 
CALL CLEAN 


^ RETURN ^ 
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d . SUBROUTINE RFR 


/sUnROUTINlA 
V -FUR J 


INITIALIZE 
BCDX, BCDY, 
BCDY1, BCDY2, 
AND BCDY3 
WRITE HEADING 
REWIND 3 


CALL CAMRAV 
CALL SMXYV 
CALL SCOUTV 
SL1 CONSTANTS 
ICT,JK AND PI 


1 

3 1 


CALCULATE 
SPPF (I), 
FIPW (I) 



COMP 

cons: 

HC,D! 
AND I 

UTE 

PANTS 

f,DX 

m 


20 } 

r 

COMPUTE 
FHZ AND 
FNW (I, J) 

49 i 


COM 

S,QX 

AND 

PUTE 
, QY 
QW 


COMPUTE HI, 
GAM2, OMEGA 

'icr - o, ic ^ i 


BEG 

LOO 

IN 21 
P 



BEGIN 40 
LOOP 



• 

20 ! 

h , 

| IA « JA H | 

J 

. 


KA 
COM 
PI PI 

JA 

PUTE 

« j 



[ PIWW -oof 

J 


BEGIN 22 
LOOP 




OMEGA - 
(F3RD (IA* 2 *Pl) 


40 

[TontinuF] 


COMPUTE 
FMl, FM2, 
RMA> KJK, 
- HMN, AC, 
EC, CC, 
POMEG, 
POWJ2, 
PIWW 



CONTINUE 

i 



CALCULATE 
OMEG (I) 
SWW (I) 

PIVVG 1 (I) 
SSSF (I) 

| SPPP (I) 



IC= 2 

FT - OMEG (l) 
G1 = SWWU) 

G2 = SSSF{I) 

G3 - PIWGl (I) 
AT SI = 00 
ATS2 = 00 
ATS 3 =00 


II 



CA LCU IA 1 1 
ATS1, AISa, 
A1S3, A I L 
A 12 Mi 
M OMI c,(f) 
GI SWW(I) 

02 * SSSI (0 

03 PfW Ol fl) | 


cm OK I OR 
MIN-M AX 
VALULS or 
SWW (I) 
SSSF(I) 
PIWGl (I) 




21 ' 


[ CONTINUE 
l KM - I 

J 


98 ^ 

1 


END FILE 3 
REWIND 3 
CALL PRNT 
SET PLOT 
SCALES 


301 


CALL SMXYV 
PLOT SPPF, 
SWW, SSSF, 
SPPP, PIWGl 
CALL CLEAN 


^ RETURN 

























SUBROUTINE RSF 


( subroutine^ 
^ rsp y 


INITIALIZE 
BCDX, BODY, 
BCDY1, BCDY2, 
BCD'* J, BCDY1 
WRITE HEADING 
UPWIND 3 


CALL CAMRAV 
CALL SCOBTV 
GALLSMXYV 
JK - 2 
1C1 -0 

PI - I 1415927 


CALCULATE 

SPPF(I) 

FIPW(I) 


COMPUTE 
HC, DX, 
EM, TCON 


20 


COMPUTE 

FNW(J,K) 


49 


COMPUTE 
S QX, QY, 
AND QW 


COMPUTE 
Hi GAM2, 
AND OMEGA 


BEGIN 21 
LOOP 



To 1 

IA = JA + l | 

U 1 

I 


KA = JA 

COMPUTE 

PIPP 


OMEGA - 


I PIWW -00 1 

T 

40 if 



TcONTINUE [ 

] BEGIN 22 1 ] 

r 


217 


COMPUTE 
XJ, XM, Y K , 
VN, FJK, 
PMN 


COMPUTE 
RMA, HJK, 
HMn, AC, BC, 
CC, POMEG, 
P0WJ2 
PIWW 



CALCULATE | 
OMEG(I) 
SWW(I) 
PlWGl(I) 
SSSF(I) 

SPPP(I) 



NO 


IC ^2 

FI = OMEG(I) 
G 1 = SOT (I) 
G2 = SSSF(I) 
G3 = PIWGI(I) 
ATS1 = 00 
ATS 2 = 0 0 
ATS3 - 0 0 


H 


41 


CALCUIAI K 
AIM ATiS2 
ATS3 All 
AT2, AT3 
FI * OM^CU) 
C,1 =SWW(I) 
C2 SSSI { l ) 
G3 PIWGI(I) 


CHICK I OH 
MIN- MAX 
VALUFS or 
SWW(I) 
SSSF(I) 
riWGi(i) 




21 J 


CONTINUf 
KM - I 


f 

9H I 

LND I 
REWlf 
CALL 
SET P 
SCALI 

I LI i 
JD I 
PR NT 
LOT 
S 


301 


CALL SMXYV 
PLOT SPPF, 
SWW, SSSF 
SPPP PIWG1 
CALL CLEAN 


^ return ^ 
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H. PROGRAM LISTING 
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0000000000004 ** 


IBFtc rand 

program random 

COMPUTER PROGRAM FOR PREDICTION OF STRUCTURAL VIBRATIONS 
DUE TO FLUCTUATING PRESSURE ENVIRONMENTS, COMBINATION 
OF PROGRAMS RSRPC1, RFRPC1, AND RSFRP1 , RECTANGULAR 
CYLINDRICAL SHELL PANEL CROSS-REINFORCED WITH RIBS AND 
STRINGERS, BOUNDARY CONDITIONS - FOUR EDGES SIMPLY-SUPPORTED, 
FOUR FDGES CLAMPED, TWO OPPOSITE EDGES SIMPLY-SUPPORTED 
WHILE OTHER CLAMPED. INPUT EXCITATION SPECTRAL DENSITY 
AND OUTPUT DISPLACEMENT, STRESS AND ACCELERATION SPECTRAL 
DENSITIES ARE PLOTTED IN GRAPHIC FORMS. 

DEVELOPED BY CHRYSLER HUNTSVILLE OPERATIONS, UNDER MSFC 
CONTRACT NAS8-21403. 

DIMENSION FNW(10flO)«F3RD(40),S3RD(40) »SPPp(400) 

DIMENSION FIPW(40I ,SPPF(40) ,OMEG(400) ,PIWGl(40n) 

DIMENSION SWW(400) ,SSSF<400) 

COMMON /INPUT/ PL,B,RHO,HS,CI ,X, Y,FINN,A1,C,PLP,BP, 

J A2,F,EP,VIP,AL1,BL1,AU, AI2»H2,HP,RAD,RH0P,RLA 
COMMON /OUTPUT/ FNW , F3RD , S3RD , F I PW , SPPF, OMEG , P I WGt , SUW , 
1SSSF,SPPP,N3R,KM, JK, ICT.PAS, ATSl,ATS2,ATS3, ATt, AT?,AT3, 
2QX,0Y,QW,GAM2, IRD 
C RFAP AND PRINT INPUT DATA 

READ (5,133) NP , N3R 
RFAD(5,5)PL»B,RH0,HS 
RE AD ( 5 » 5 ) C I ,X,Y,FINN 
READ (5,5) A1 , C , PIP , BP 
READ (5,5) A2 » E , EP , V I P 
READ (5,5) AL1.BL1* AH, A 12 
READ (5,5) H2, HP , R AD , RHOP 
READ (5,8) (F3RD( I ) , S3RD < I ) , I»1,N3R) 

WRITE (6,103) PL,8,RH0, HS,CI , X , Y , F I NN , A1 , C 

WRITE (6,132) PLP,8P, A2,E,EP, VIP,AL1,BL) , AH, AI2,H2,HP,RAD,RH0P 
GO TO <200, 201, 202, 200), NP 

C CALCULATING RESPONSFS OF FOUR EDGES SIMPLY-SUPPORTED PANEL 

200 CALL RSR 

IF (NP .NE, 4) GO TO 99 

C CALCULATING RESPONSES OF FOUR EDGES CLAMPED PANEL 

201 CALL RFR 

IF (NP .NE, 4) GO TO 99 

C CALCULATING RESPONSES OF TWO OPPOSITE EDGES SIMPLY-SUPPORTED 
C WHItE OTHER TWO CLAMPED PANEL 

202 CALL RSF 
99 STOP 

5 FORMAT ( 4F15 , 8 ) 

8 FORMAT (8F8.1) 

103 F0RMAT(1H1,10X,16HINPUT DATA ,//, 

120H PL = PANEL LENGTH =,E12,5,/, 

218H B a PANEL WIDTH a,El2.5,/, 

330H RHO = MASS DENSITY OF PANEL a,E12.5,/, 

423H HS = PANEL THICKNESS »,E12,5,/, 

521H Cl = DAMPING RATIO * t E%2,5,/, 

629H X s COORDINATE OF VECTOR R s,El2.5,/, 

729H Y s COORDINATE OF VECTOR R a,E!2,5,/, 

862H FINN a ONE NTH OCTAVE FREQUENCY INCREMENT = 

9,E12,5,/,34H A1 a DECAYING CONSTANT - LENGTH a E12.5,/, 
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*5 2 1 H C s SPEED OF SOUND = E12.5) 

132 FORMAT ( 48H PLP = LENGTH OF PANEL SUBJECTED TO EXCITATION =612,5,/ 
146H BP = WIDTH OF PANFL SUBJECTED TO EXCITATION =E12.5,/, 

233H A2 = DECAYING CONSTANT - WIDTH sEl2.5>/. 

330H E = YOUNGS MODULUS OF PANEL =E12.5,/, 

436H EP = YOUNGS MODULUS OF STIFFENERS =£12.5,/, 

523H VIP = POISSONS RATIO =E12,5,/, 

644H ALL * SPACING OP Y - DIRECTION STIFFENERS = E]2.5,/, 

744H BL1 s SPACING OF X - DIRECTION STIFFENERS =Fl2,5,/, 

857H All s MOMENT OF INERTIA OF ONE Y - DIRECTION STIFFFNER =£12,5 
9/.57H A 1 2 s MOMENT OF INERTIA OF ONF X - DIRECTION STIFFENER = 
$E12.5,/,43H H2 = LARGEST HEIGHT OF STIFFENERS AT VECTOR R =F12,5, 
$/,43H HP s SMEARED-OUT THICKNESS OF STIFFENERS =E12.5,/, 

S24H RAD s RADIUS OF SHELL =E12,5,/, 

T36H RHOP = MASS DENSITY OF STIFFENERS =E12,5,///) 

133 F ORMAT (415) 

END 

sibftc pnt 

SUBROUTINE PRNT 

DIMENSION FNW(IO.IO) , F3RD < 40 ) , S3RD ( 40 ) * SRPP ( 400 ) 

DIMENSION FIPW(4Q),SPPF(40>»0MEG(400) ,PIWG1 (400) 

DIMENSION SWWM00 5 *SSSF( 400 5 

COMMON /INPUT/ PL,B,RHO,HS,CI ,X, Y»F!NN,A1,C,PLP,BP. 

1A2.E.EP, VIP. AL1.BL1, All, AI2.H2.HP, RAD.RHOP.PLA 
COMMON /OUTPUT/ FNW , F3RD, S3RD , F I PW , SPPF , OMEG , P I WQ1 , SWW, 

1SSSF ,SPPP, N3R.KM, JK ,IGT,PAS»ATSi » ATS2 , ATS3 » ATI , AT2 » AT3, 
2QX.QY.QW.SAM2, I RD 
PIs3, 1415927 
WRITE (6,110) 

DO 200 1=1, N3R 

200 WRITE (6,115) F3RD ( 1 ) , S3RD ( I ) , SPPF ( I) , F I PW ( I ) 

WRITE (6,100) 

DO ?01 J=1 , 9 
DO 201 K=1 , 9 
FHZ s FNW( J,K)/<2.*PI ) 

201 WRITE (6,120) J.K.FHZ ,FNW(J,K) 

WRITE (6,125) 

206 WRITE (6,130) ( OMEGU ) . SWW ( I ) , SSSF ( I) , P I WG1 ( I) , SPPR ( I ) , I *1 , KM > 

207 IF (JK ♦ EQ , 2) GO TO 203 

WRITE (6,102) ATS1, ATS2,ATS3 
WRITE (6,101) ATS , AT2,AT3 

WRITE (6,131) QX,QY.QW,GAM2>PINN 
GO TO 204 

203 WRITE (6,102) ATS1.ATS2.ATS3, PAS 
PAsSQRT(PAS) 

WRITE (6,101) ATI » AT2, AT3, PA 

WRITE (6,131) QX,QY,QW,GAM2,FINN 

204 WRITE (6,105) 

IF (ICT ,EQ, 0) GO TO 208 
DO 202 1=1, ICT 

READ (3) POMEG , J, K , M , N , P0W J2 

202 WRITE' (6,104) POMEG , J , K , M , N , POW J2 
REWIND 3 

203 WRITE (6,50) 

RETURN 
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50 Format ( i2hiend of data) 

100 FORMAT ! lHl , 36H NATURAL FREQUENCIES //.2X4MJ K,l 

10X3HFMZ,20X3HFNW/1 6X2HHZ , 18X7HRAD/SEC// ) 

101 FORMAT ( 13H RMS VALUE ,4E16.8,/) 

10? FORMAT!/, 13H MS VALUE ,4E16.8,/> 

104 F0RMAT(F12,5,4I3,3E55,3) 

105 Format ( iihi frequency, 5x,4hmode,5x,i6hjoint acceptance,/, 

14x,5NHERTZ,6X,7HINDICES,8X,6HSQUAP6./,14X,10HJ K M N , / ) 

107 FORMAT<11X,30HOVERALL MEAN-SQUARE PRESSURE =,Ei5,7) 

110 FORMAT!///, 38X16HINPUT EXC I TAT I ON//AX14HONE-3RO OCTAVE. 7X?0MOVfc-3R 
ID OCTAVE LEVEL, 7X16HSPECTRAL DENS I T Y , 7X1 6HSP6CTR AL DENS I TV/flXl 4WME 

2 AN FREQUENCY, 12XHHIN DEC I BELS , 11X17H I N DEC I BELS/HgRTZ , 6X1 7H I N »Si 

3 SQ/RAD/SEC/14X4HF3PD,19X4HS3RD.22X4HSPPF,20X4HFIPW//) 

115 FORMAT(F20,3,3XF20.3,6XF20.3,6XE20 .5) 

120 FORM AT ! 2 1 3 , Ft5.5,7X,Fl5.5) 

1?5 F ORMAT ( 13H1 FREQUENCY , 3X , 1 2ND I SPLACEMENT , 7X , 6HSTRESS , 7X , 
11?HACC6LERATI0N,5X, JOHEXCITATION, /,18X,8HRESPONSE, 4X, 

216H RESPONSE , 4X » 8HRESP0NSE » 

3/,6X,5HHERTZ,3X.16H INCH SQ/HERTZ ,2X,12HPSI SO/HERT7 , 5X , 10HG SQ/ 
4HFRTZ,4X,14HPSI SQ/RAD/SEC , / ) 

I3n FORMAT! 1X,F12.5.4E16,8) 

131 FORMAT!/, 5H QX =E12 , 5 , 2X * 4HQ Y aF12 , 5 , 2X , 4HQW a£12 . 5 « 2X , 8HGAMMA2 s 
1E12.5,/,8H FINN =,F5,l/> 

END 

SIBFTC GRID 

SUBROUTINE GRID1V (L , XL , XU , Yt , YU , DX , DY , NN , MM , H , J J , NX, NY ) 

C GENERAL GRID SUBROUTINE FOR SC4020 — LINEAR OR NONLINEAR IN - 

C EITHER VERTICAL OR HORIZONTAL 

DIMENSION XYLf 2) ,XYU!2) ,DXY(2) , NM < 2 > , I J(2> , NXV <? ) , XYLAB<2> . UR£FR< 
12) ,K ! 2 ) » LLL ! 2 ) ,LUU(?) ,LL(2) . LU ( 2 ) , MU < 2 ) , MTU ! 2 ) , ML ! 2 ) , MU ! 2 ) , I TT ! 2 ) 
2 , I SPACE ! 2 ) , IT0P!2) 

EQUIVALENCE ! I TOP < 1 ) , MUH ) , ( I TOP ( 2 ) , MUV J 
C STOP TYPE--MAY NOT MEED 

NYs-3 

70 CALL STOP-TV 
CALL BRITEV 
IF(|.-2) 95.110,95 
95 INDsL-1 

100 CALL FRAMEVMND) 

110 CALL SETMOV ! MTU l ) ,MTU ( 1 ) , MTL ( 2 ) , MTU ! 2 )) 

CALL SETCOV (IWIDE, IHJGH) 

CALL MSXYV ! K ( 1 ) , K ! 2 ) ) 

120 XYUDsXL 
XYL(2)=YL 
XYU < 1 ) *XU 
X YU ( 2 ) s YU 
DXY<1)«DX 
DXY ( 2 ) sDY 
NM(1)=NN 
NM(2)sMM 
IJ( 11*11 
I J < ? > = J J 
NXY(1)=NX 
NXY(2)sNY 
DO 230 Is 1,2 
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ISPACF ( I ) = 0 
LUU ( I ) s 0 

IF(NXY(I>) 229,230,230 

229 NXY(I >=7-NXY( I ) 

230 CONTINUE 

IF(NXY(2) >240,244,240 
24 n ISPACE(1 )sNXY<2)*IWIDE+IW!DE*6 
UUU(2)sIHIGH+6 

244 [F(NX> 245,250,246 

245 ISPACF(2)*IHIGH/2 

246 ISPACE(2)sISPACE{2)+IHIGH+6 
LUU(l)«(NXY<l)*IMlDF)/2 

250 CALL HOLDOVUND) 

c compute label margins, if ind is zero 

30D DO 445 1=1,2 

LLL < I ) = ISPACEt I) 

XYLA8( I )=XYL< I ) 

IF ( - I J ( I ) ) 330,360,360 
330 IF(XYL(I)) 340,360,340 
34n IF(XYL< I)*XYU<I>) 350,350,360 
350 XYL AB ( I > = 0 .0 
355 LLL( I )=0 

360 IF (NXY(l) ) 370,380,370 

370 LLL(i>sMAX0ailUtl) + TWIDE + 2»LLL<l> > 

380 IF (IND) 420, 3.90,420 
39o LL(I)*LLL<n 
' LU(I)cLUUU) 

420 ML ( I ) = LL( I ) *MTL f I ) 

MLIU )=LU( !)*MTUM > 

440 I TT ( I )=MLM )+MU( I > 

445 LREFR ( I ) sMTL ( I ) +3 
C END OF LOOP FROM 300 

460 IF (MINO(NMIl) ,NM(2) ) >465,660,660 
465 ITsHAXOMTTU) , ITT<?) ) 

480 MU(1) = IT-MLU> 

MH(2)=IT-ML(2> 

C ERROR TESTS AND SCALE 

660 IF(K(1) >662,669,662 

667 CALL ERRNLV(XYL(l)rXYU(l),ML(l),MU(j),DXY(l>) 
XVLA8{l)sXYL{l> 

GO TO 670 

669 CALL FRRLNV< XYL(l) , XYU<1 ) ,ML( 1 > ,MU{ 1) , DXY(l) > 

670 CALL XSCALV ( XYL ( 1 ) , XYU ( 1 ) , ML ( 1 > , MU { 1 ) ) 

IF ( K ( 2 > ) 682,692,682 

687 CALI ERRNLV ( XYL ( 2 ) , XYU < 2 ) , ML ( 2 ) , MU ( 2 ) , DXY ( 2 ) > 
XYLAB(2)=XYL(2> 

GO TO 700 

69? CALL ERRLNV ( XYL ( 2 ) , XYU < 2 > , ML ( 2 ) , MU ( 2 > , RXY ( 2 ) > 
700 CALL YSCALV ( XYL < 2 ) , XYU ( 2 > . ML ( 2 ) , MU ( 2 > ) 
ITT(15=NXV(XYLAB(1) >-!SPACE(l> 

ITTf 2)=NYV(XYLAB(2> >-ISPACE(2) 

DO 790 1=1,2 

I TOP ( I >*1023-MU< I > 

758 IF( ITT< I )-MTL( I) >760,780,780 
760 I F f - 1 J f I ) > 770,790,790 


4a 



oooorioooooooooooooo 


77n U(i)s-iJU) 

GO TO 300 

780 LREFR ( I ) = I TT { I ) +2 
790 CONTINUE 

BOB LREFR(i)=LREFR{i)+IWIDE+IWIDE/2 
LREFR(2)=l.REFR{2) + IHIGH/2 
810 I F { K ( 1 ) > 820.850,820 

820 CALL N0NLNV<1,LREFR<2),ML(2).MUV,XYL(1>,XYU(1WDXV<1 ), NMC1),IJ<1> 
l, NX, I W I DE ) 

GO TO 870 

850 CALL L]NRV(1,LRFFR{?> ,ML(2) . MUV . X YL < 1 > , XYU <1 ) , DXY ( 1 ) , NM ( 1 ) , I J < 1 ) , 

INX. IWIDE) 

87n I F ( X { 2 ) > 880,910,880 

880 CALL N0NLNV<2,LREFR(1) ,ML(1) , MUH , XYL ( 2 ) , XYU< 2 ) , DxY ( 2 ) , NM < ? > , I J < 2 ) 
l.NY, IHIGH) 

900 GO TO 924 

910 CALL LINRV(2,LREFR(1> ,ML(1) , MUH.XVLI2) , X YU < 2 ) , QX Y ( 2 > , NM < 2 > , I J 1 2 ) , 

INY, IHIGH) 

C TO CLEAR ERROR INDICATOR —IF GRID DATA IN ERROR .ERROR INDICATORS 

C 1MIGHT NOT BE CLEARED 
924 K(1)=NXV(XYLI1) > 

K(2)=NYV(XYL(2) ) 

930 RETURN 
END 

SORIGIN A 

$ I BFTC RSR1 

SUBROUTINE RSR 
PROGRAM NUMBER 823-1002-6 
RSRPCl 

CALCULATION OF RESPONSE of SIMPLY SUPPORTED RECTANGULAR PANELS 
CROSS-REINFORCED WITH ST I FFENFRS 
INPUT PARAMETERS 
PL = PANEL LENGTH 
B s PANEL WIDTH 
RHO = MASS DENSITY OF PANEL 
HS = PANEL THICKNESS 
Cl = DAMPING RATIO 
X = COORDINATE OF VECTOR R 
Y = COORDINATE OF VECTOR R 
FINN = ONE NTH OCTAVE FREQUENCY INCREMENT 
A = CONSTANT 
C = SPEED OF SOUND 

PLP = LENGTH OF PANFL UNDER EXCITATION 
BP = WIDTH OF PANEL UNDER EXCITATION 
F3RD s ONE-THIRD OCTAVE BAND CENTER FREQUENCY 
S3RD s SOUND PRESSURE LEVEL 

DIMENSION FNWem.lO) , F3RD < 4ft ) , S3RD < 40 ) ,BCDx<12) , BODY (12) 

DIMENSION F I PW ( 40 ) , SPPF ( 40 ) , OMEG ( 400 ) , P I WGl < 40 0 ) ,RCDYK1?> 
DIMENSION BCDY2(12),SWW(400) ,SSSF( 4 nO) ,RCDY3(12) 

DIMENSION SPPP ( 40Q ) 

COMMON /OUTPUT/ FNW , F3RD , S3RD , F IPW , SPPF , OMEG , P J WGl , SWW , 
1SSSF,SPPP,N3R,KM, J«, ICT,PAS,ATS1,ATS2,ATS3,AT1,AT?,AT3, 
2QX,0Y,QW,GAM2, I RD 

COMMON /INPUT/ PL , B , RHO , HS , C I , X , Y , F I NN , A1 , C , PLP , BP , 
1A2,E,EP,VIP,AL1,BL1,AI1,AI2,H2,HP,RAD,RH0P,RLA 
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DATA RCDX/72HFREQUE W CY <H?> 

1 / 

DATA BCDY/72HACCEIERATI0N G SQ/H* 

5 / 

DATA BCDY2/72HDISPLACEMENT INCH SQ/HERTZ 
1 / 

DATA BCDY3/72MSTRESS SPECTRAL DENSITY PSI SQ/HERTZ 
I / 

DATA BCDY1/72HSPPP DB/HZ 
1 / 

WRITE <6,1083 
!CK = 1 

1001 REWIND 3 
JKsl 

Pls3, 1415927 
SWMI = 1.E10 
SWMX= X.E-10 
SSNIs 1.E10 
SSHXs l.E-10 
PIMIs 1.E10 
PIMXs l.E-10 
CALL CAMRA V ( 9 ) 

CALL SMXYVd.O) 

C CALCULATE EXCITATION SPECTRAL DENSITY 

DO 9 Isl,N3R 

SPPPU)sS3RD( n-l0,»AL0G10(0 r 2315»F3RD( I ) ) 

FIPW{ nsl t /«2.0*PI)«J.0,O**(CSPPr< I)-170.976)/10.0) 

9 continue 

C CALCULATE natural frequencies 

HCsF*H$*HS*HS/(12.*<1,-VIP*VIP) ) 

DY=HC+EP*AI2/AL1 

DX=HC+EP*AI1/BL1 

BM = RHO*HS+RHOP*HP 

DO 20 I «1 » 9 

DO 20 Jsl » 9 

RMs I 

RNs J 

FHZ = DXMRM /PL)*#4 + 2.*HC#<RM*RN/(PL*BM**2 + DY«<RN/8)**4 

1+E*HS/(RAD*RAD»(PI )#*4*<1 ,*<PL«RN/(B*RM) }#*2)*#2) 

FNWU , J)=PI*PI*SQRT( 1 , /8M ) »SQRT { FHZ ) 

20 CONTINUE 

C CALCULATE CONSTANT TO CHANGE DISPLACEMENT INTO STRESS 

QXsO . 0 
QY = 0 , 0 
QWsO , 0 
S=BP*PLP 
DO 49 I J»1 » 5 , 2 
DO 49 IKsl,5.2 
RMs I J 
RN=IK 

FSINsSIN(RM*P!*X/PL)*SIN(RN*PI#Y/R) 

QC0NsRM«RN*<'DX*(RM/PL)**4*2.*HC*<RM*RN/<PL*Bn*»2 + DY#<RN/En**4> 
QX=OX+( (RM*PI/PL)«*?+V!P#<RN#PI/B3**^)*(FSIN/0C0N) 

QY=QY+< <RN*P1/B)**2*VIP*(RM*PI/PL>**2>*<FSIN/QC0N> 

49 QWsQWt-FSlN/QCON 
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HisHS+H? 

GA.M2sE»H1*H1*E/<4 ,*(!.- VIP* VIP )**2>*< ( QX*OX+QY#Oy ) / <QW*0W > ) 

C CALCULATE EXCITATION PRESSURE PROM ONE-THIRD OCTAVE DATA 

inOO CONe 2. 0*»(1 ,0/FINN) 

OMEGA=F3RD(1)*2.0*PI 

I CT = 0 

ICsl 

DO 21 1=1,400 
N3R1 =N3R-1 
DO 40 JAal , N3R1 

IF ( OMEGA , GT , ( F3PD < N3R ) *2 , *P I ) ) GO TO 99 

IF ( OMFGA -< F3RD ( JA)*2,*PI ) >25,25,30 
25 KAaJA 

GO TO 35 
30 IAsJA+1 

IF (0MEGA-<F3RD( I A ) *2 , *P I > >25,25,40 

35 PIPPsFIPW(KA)*(FIPW(KA+l)-FIPW(KA> ) * ( OMFQA-F3RD ( KA > #2 . *P I ) / ( (F3RD 
l(KA*l)-F3RD(KAn*2,*PI ) 

GO TO 45 
40 CONTINUE 

45 PIWWsO.O 

DO 22 J=1 » 5 
DO 22 K = 1 . 5 
DO 22 M = 1 * 5 
DO 22 N=1 , 5 
RK = K 
RNsN 
RJ = J 
RMsM 

FMl=SIN<RJ*PI*X/PL)*StN(RK#PJ*Y/B) 

FM2=SIN(RM#PI»X/PL)*SIN(RN*PI*Y/B) 

RMA=RHO*HS*B*PL/4.0 

HJK=ABS<1,/<RMA*SQRT( {FNWC J»K)#FNW< J,K)-OMEGA«OMEGA>**2 
1*(2.*CI*FNW( J,K)*OMEGA)**2) ) ) 

HHNbABSU. /(RMA*SQRT< (FNWIM, N)*FNW(M,N)-0MEGA*0MEGA)**2 
l+(2.»CI»FNW(H,N)*OMFGA)»*2) ) ) 

C CALCULATE JOINT ACCEPTANCES 

AC=a.-(OMEGA/FNW( J,K> )*#2)»(1,-(0MSGA/FNW(M,N))#*2)+4,*CI»CI^0MEG 
1A**2/<FNW( J.KUFNW(M,N) ) 

BC = -2.*( (CI*OMEGA/FNW( J,K))*U ,-*(OMEGA/FNW<M,N))**2 I -Cl* OMEGA/ 
1FNW<M,N)»(1.-(0MEGA/FNW( J,K> )**2> ) 

CCsAC*AC*BC*BC 

CONl5(Al»OMEGA#PL/C)*(Al*OMEGA*PL/C) 

C0N2*R J*P I 
CON3=RM»PI 

C0N6 = C0N2*C0N3*( ( 2 . +EXP < - Al«PL*OMEGA/C ) # { < -1 . ) »« ( J+l > ♦ ( ( -1 , ) #* 
i(M+l> ) ))/( (CON1+CON2#CON2)*(CON1+CON3*CON3) ) J 
CON7*RK*P I 
C0N9sRN»PI 

IF (J ,NE. M) GO TO 46 

PJJM = (Al#PL»OMEGA/C/2.*(l ,/(C0Nl*C0N2*C0N2)+ 
ll,/<CONl + CON3*CON3n+CON6> 

GO TO 47 

46 PJJM a C0N6*RM / ( CON1+CON3#CON3 >* ( ( ( -1 , > ** I J-M)-l . ) / ( 2 . * 
l(RJ-RM>) + ( <~1, )*«< J*M>-1, >/<2.*(RJ*RM> > ) 
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2 + RJ /(C0Nl + C0N2*CnN2)*( ( (-1. . )/< 2 , »<RM-RJ> ) 

3+< (-1. W«2.*(RM+BJ) ) > 

47 C0N1*(A2*0MEGA*F /C ) * ( A2*0MFGA*R /C) 

CON6sCOM7*CON8*< (2,4EXP(-A2#B *OMEQA/C)*( t-1 . ) **<K*1 > + < <-t . > 

1 (N + l) > } )/( (CONI +CON7#COW7 )*( CONI +C0N8*C0N8 ) ) ) 

I F ( K .ME. N> GO TO 48 

PJKN = ( A2*PL«nMEGA/C/2.Ml t /(CONl + CON7*CON7> + 

1 1 . / (CON1+CON8*CON8 )> +C0N6 } 

GO TO 23 

4H PJKN a CQN6*RN /(CON1+CON8»CON8)*< ( <'l, )/<2.* 
l(RK-RN) )+<<-!. >**<K*N)-i, )/ <2 . * ( RK + RN ) ) ) 

2 +RK /(C0N1+C0N7*C0N7)«( «(-l. )/{2,*(RN-RK») 
3+U-i. )**(N*K>-1, )/<2,*<RN+R«) ) ) 

23 P0MEG=0MEGA/(2,«PI ) 

P0WJ2=PJJM*PJKN*AC/SQRT<CC) 

PIWWsPIWW4FM1*FM2*HJK#NMN#P0WJ2 
IF (PONEG , GT . 5, ) GO TO 22 
WRITE (3) POMEG , J , K , M , N , POW J2 

ICT=ICT+1 

22 CONTINUE 

OMEG ( I ) = OMEGA/(2,*PI ) 

CALCULATE DISPLACEMENT, STRESS, AND ACCELERATION 

PIWWaPIWW*S*S*PIPP 

SWWU )*PIWW*2.*PI 

PIWG =OMEGA**4*PIVW 

PIWGK I JaA.SlSORSE-OS^PIWG 

PSSW=GAM2«PIWW 

SSSF( I >=2.*PI*PSSW 

SPPPU 5 = 2. *Pl*P IPP 

CALCULATE MEAN SQUARE 

IF ( IC ,NE. 1 ) GO TO 41 

FlaOMEG ( I ) 

GlaSWWU ) 

G2aSSSFU ) 

G3aPlWGl(I) 

ATSlaO ,0 
ATS2aO ,0 
ATS3s0 t 0 
IC = 2 

GO TO 21 

41 ATS1=ATS1+(G1*SWW( I ) ) /2 . * ( OMEG ( I ) ~F1 ) 

ATS2aATS2*(G2+SSSF( I > ) /2 , * { OMEG { I ) -Fl ) 

ATS3aATS3+(G3 + PIWGl< I )) /2 . * <OMEG( M-Fl ) 

CALCULATE ROOT-MEAN SQUARE 
FlaOMEGM ) 

GlaSWWU ) 

G2aSSSF( I ) 

G3 = PIWG1U ) 

IF(SWWU) , LT . SWM I ) SWMIaSWW(I) 

IF(SWWU) , GT ♦ SWMX) SWMX = SWWU) 

IF(SSSFU) , LT , SSMT)SSMlaSSSF(I) 

IFISSSFU) .GT, SSMXISSMXsSSSFU ) 

IF(PIWGKI) ,LT. P I M I ) PlMIapiWGiU) 

IF ( P I WG1 < I > .GT. PIMX) PIMXaPIHGIM) 

OMEGAa OMEGA* CON 
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21 CONTINUE 
KM* I 

GO TO 98 
9Q KMs I -1 
9fi END FILE 3 
REWIND 3 
ATlaSGRT ( ATS1 ) 

AT2pS0RT(ATS2) 

AT3=S0RT( ATS3) 

CALL PRNT 
DSW a SWMX/1.E9 
DSS a SSMX/1.F9 
DPI a PIMX/1.E9 
DO 301 Ial.KH 

IF(SWW(I> , LT . DSW ) SWW(I)aDSW 
IF(SSSFd) .LT, DSS) SSSFd)«DSS 
IF(PIWGKI) ,LT. DPI) P I WG1 (I ) aDP I 
301 CONTINUE 

C PLOT EXCITATION, DISPLACEMENT, STRESS, AND ACCELERATION 

CALL 0UIK3VI-1 ,44,8CDX, BCDY1, -N3R , F3RD , SPPF ) 

CALL SCOUTV 
CALL SMX YV ( 1 » 1 ) 

CALL OUIK3V(-1,44,BCDX,8CDY2,-KM,OMEG,SWW> 

WRITE (16,106) ATI , F I NN 

CALL QUIK3V<-l,44,BCDX,BCDY3»-KM t 0MEG,SSSF ) 

WRITE (16,106) AT2 » F I NN 

CALL 0UIK3V(-1,44,BCDX,8CDY,-KM,0M6G,PIWG1 ) 

WRITE (16,106) AT3 ,F J NN 
CALL CLEAN 

106 FORMAT (11X,31H ROOT-ME AN-SOUARE RESPONSE a E15 , 7 , 5X , 6Hp I NN = 
1F6,2,10X,6HRSRPC1 ) 

108 FORMAT (96H1CALCULAT ION OF RESPONSE OF SIMPLY SUPPORTED RECTANGULAR 
1 PANELS CROSS-REINFORCED WITH STIFFENERS) 

RETURN 

END 

SORIGIN A 

SIBFTC RFR1 

SUBROUTINE RFR 
PROGRAM NUMBER 823-1002-7 
RFRPCl 

DYNAMIC RESPONSE OF FOUR-SIDE FIXED RECTANGULAR 
SHELL PANELS CROSS-REINFORCED WITH STIFFENERS 
INPUT PARAMETERS 
PL * PANEL LENGTH 
B a PANEL WIDTH 
RHO a MASS DENSITY OF PANEL 
HS a PANEL THICKNESS 
Cl a DAMPING RATIO 
X a COORDINATE OF VECTOR R 
Y a COORDINATE OF VECTOR R 

FINN a CONSTANT TO DETERMINE OCTAVE BAND CENTER FREQUENCIES 
A1 e CONSTANT 
C a SPEED OP SOUND 

BP a LENGTH OF PANEL UNDER EXCITATION 
PLP a WIDTH OF PANEL UNDER EXCITATION 
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C F3RD = ONE-THIRD OCTAVE BAND CENTER FREQUENCY 
C S3RD s SOUND PRESSURE LEVEL, 

DIMENSION FNWdO, 10 > , F3RDM0 ) , S3RD(40 > » RCDX( 12) »BCDY(12 > 
DIMENSION FIPW(40),SPPF(40)> 0MEG(400) ,PIWGi(400>,BCnYl<J2> 
DIMENSION BCD Y2 ( 12 ) , SWW ( 400 > , SSSF < 400 > ,BCDy3<1 2 ) ,BCDY4(J?) 
DIMENSION SPRP (400) 

COMMON /OUTPUT/ FNW , F3RD, S3RD. FIPU, SPPF, OMEG,P f WGl , SWW, 

J- SSSF ,SPPP»N3R»KM»JK> ICT,PAS»ATS1,ATS2,ATS3,AT1.AT2,AT3, 
2QX,QY,QW,GAM2, IRD 

COMMON /INPUT/ PL , B , RHO, HS, C J , X , V , F INM, A1 , C, PLP, BP, 
1A2,F,EP, VIP, AL1,BL1,AU, AI2,H2,HP,RAD,RH0P,RLA 
DATA BCDX/72HFREQUENCT (HZ) 

1 / 

DATA BCDY/72HACCELERATION G SQ/HZ 
1 / 

DATA BCDY2/72HDISPLACEMENT INCH SQ/HFRTZ 
1 / 

DATA BCDY3/72HSTRESS SPECTRAL DENSITY PS I SQ/HERTZ 
1 / 

DATA BCDY1/72HSPPF DB/HZ 
1 / 

DATA 8CDY4/72HSPECTRAL DENSITY PSI SO/HERTZ 
1 / 

WRITE (6,108) 

REWIND 3 

JKfi? 

I C T a 0 

PI=3, 1415927 
SWM I = 1 , E10 
SWMX= l.E-10 
SSM I = l.EtO 
SSMXs l.E-10 
PIMI* 1.E10 
P1MX= l.E-10 
CALL C AMRA V ( 9 ) 

CALL SMXYV (1,0) 

C EQUATIONS 11 AND 12 

C CALCULATE EXCITATION SPECTRAL DENSITY 

DO 9 !sl,N3R 

SPPF ( I )sS3RD( I )-10,*ALQGl0(0,2315*F3RD( I ) ) 

FIPW( I )=1,/(2,0*PI ) *10 , ( (SPPF ( I >-170.5765/10, 0 ) 

CONTINUE 

CALCULATE NATURAL FREQUENCIES 
EQUATIONS 1 AND 2 
HC=E*HS*HS*HS/(12,*(1,-VIP*VIP) ) 

DYsHC+EP*AI2/ALl 
DXsHC + EP«AU/8Ll 
BM a RHO*HS+RHOP*HP 
FCON = PI *P I /SORT ( 8M ) 

DO 20 J=1 , 9 
DO 20 K=1 , 9 
RJ u J 
RK r K 

IF(J .NE. 1 .AND. K .NE. 1) GO TO 18 
IF CJ .EG, 1 .AND, K ,N6, 1) GO TO 16 
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IF(J ,NE, \ , AND , K .Fq, 1) Go TO 17 

FHZ = DX* < 1. 5056/PL > **A + DY# ( i . 5056/R ) *»4+? . »HC» f 1,24 6 6/ <PL*»> >•* 2 
1+E*HS/(RAD#RAD*(PI > **4* ( l , + t PL /R ) * ( PL/8 ) ) « *2 ) 

GO TO 19 

16 FHZ = DX*(i ,5056/PL)*M + DY*< < RK+ . 5 ) /S ) *»4 
1+2,#HC*< <1 ,2466*(RK*,5>*( < RK+ . 5 ) - . 6366 ) 5 / < PL*PL*B*B > > 
2+E*HS/£RAD*RAD*<P! ) »*4# ( 1 , + ( ( RK+ . 5 ) /1 .50 56 >**2* ( PL/P ) * (PL/B) )«*? ) 

GO TO 19 

17 FHZ = DX*(<RJ+.5)/PL >**4+DY*<l, 5056/8 )#M 
1+3.*HC#(1,2466*CRJ*.5)*( t R J+ , 5 ) - , 636 6 ) / ( PL*PL*8*B ) ) 
2*E*HS/(RAD*RAD*(PI )#*4*(i t +(l,5056/(Rj+.5> )**2 *<Pl/B)*(PL/B>)**2> 

GO TO 19 

18 FHZ = DX*( (RJ«-.5J/PL)#*4 + DY*< ( RK+ , 5 > /R > **4 

1 + (RJ* ,5>*<RK*.5J*C «RJ*.5)-.6366>*( (RK+.5>-.6366)/<PL*PL*B*B 

?) )+E»HS/(RAn*RAD*(PI ) #* 4 *<!.+( (RK* ,5)/<RJ+.5> )*«2 
3*<PL/B)*<PL/B> >*«2> 

19 FNW(J,K> = FCON»SORT(FHZ) 

?n CONTINUE 

C CALCULATE CONSTANT TO CHANGE DISPLACEMENT INTO STRESS 

C EQUATION 16 

QXsO *0 
QYsO . 0 
OWsO.O 
SsBP*PLP 
DO 49 IJsl.5,2 
DO 49 I K = 1 1 5 > 2 
RMs I J 
RN= I K 

FSINsS{n(RM*PI*X/PL)*SINC RN#P I»Y/B) 

QC0NsRM*RN*CDX«(RM/PL>** 4*2,*HC*(RM*RN/<PL*B> >*»2+DY*(RN/p)»*4) 
GX=OX+£ (RM#PI/PL)*#2+YIP«<RN*PI/B>**2)*<FSIN/QC0N) 

OYsQY+( <RN»PI/B>**2*VIP*(RM#Pr/PL5**2)«(FSlN/OCON) 

49 0WsOW+FSIN /QCON 
Hl=HS+H2 

GAMa»E*Hl#Hi # E/<4.*<l»-VlP#VJP)**2>*< (QX*QX+GY*QY>/(QW»QW) > 

C CALCULATE EXCITATION PRESSURE FROM ONE-THIRD OCTAVE DATA 

1000 CON* 2. 0**fi. 0/FINN) 

OMEGA=F3RD(1)»2.0*PI 

ic«i 

DO 21 1=1,400 
N3R1=N3R-1 
DO 40 JA=1,N3R1 

IF < OMEGA ,GT, ( F3RD ( N3R > #2 , *P I ) > GO TO 99 
IF(0MEGA-(F3RD(JA)#?,#PI) >25,25,30 
25 KA=JA 

GO TO 35 
30 I As JA + 1 

I F(OMEG A- ( F3RD( IA)*?.*PI >>25,25,40 

35 PIPP=F[PW(KA)*(FJPW(KA+1 >-FlPW(KA) >#<0MEGA-F3RD<KA)*2,#PI >/( (F3RD 
1 ( KA*1 ) **F3RD ( KA ) > *2 , #PI > 

GO TO 45 
40 CONTINUE 

45 PIWWso.O 
DO 22 Jcl , 5 
DO 22 K = 1 ,5 


47 



DO 22 M =1 , 5 

DO ?2 N = l,5 

RJsJ 

RKsK 

RMsM 

RNrN 

C CALCULATE THE NORMAL MODE 

C EQUATION 4 

CXI = 1.5056*PI*X/PL 

CY 1 s X .5056»PI»y/9 

CXJ s (RJ+,5)*PI*X/PL 

CXM s (RM+,5)#PI*X/PL 

CYK s ( RK+ , 5 ) *P I * Y/R 

CYN 3 (RN+,5)*Pl*Y/e 

CHX1 = (EXP(CX1)+EXP(-CX1> )/2. 

SHXJ = f EXP(CXl>-EXP(-CXin/ 2 . 

CHY1 = (EXP(CY1)*EXP(-CY1) )/2. 

SHY1 3 <EXP(CY1)-EXP(-CY1) )/2. 

CHXJ = <EXP<CXJ1+EXP(-CXJ) >/2. 

SHXJ = (EXP{CXJ)-EXP(-CXJ) )/2. 

CHXM = (EXP(CXM)+EXP(«-CXM))/2. 

SHXM s (EXP(CXM)-EXP(-CXM) >/2. 

CHYK = (EXP(CYK)+EXP{-CYK) )/2. 

SHYK s (EXP<CYK)"EXP(-CYK) >/2. 

CHYN 3 (EXP(CYN)+EXP(-CYN) )/2. 

SHYN s (EXP(CYN)-EXP(-CYN) )/2. 

IF<J .GT. 1) GO TO 210 

XJ s CHX1-C0S(CX1)-,9825*(SHX1-SIN(CX1 > ) 

GO TO 211 

210 XJ s CHXJ - COS(CXJ)-(SHXJ-SJM(CXJn 

211 IF <K , GT i 1) GO TO 212 

YK = CHY1-C0S(CY1)-,9825*(SHY1-SIN(CY1) ) 

GO TO 213 

212 YK 3 CHYK«COS(CYK)-ISHYK-,SIN(CYKn 

213 I F ( M .GT, 1 ) GO TO 214 

XM s CHX1-C0S(CX1>-,9825»<SHX1-SJN<CX1> > 

GO TO 215 

214 XM s CHXM-rCOS<CXM)'(SHXM-SIN{CXM) ) 

215 I F { N , GT , 11 GO TO 216 

YN s CHYl-C0S(CYl)-.9825«(SHYl”SlN(Cyin 
GO TO 217 

216 YN s CHYN-COS<CYN)-<SHYN-SIN<CYN) ) 

217 FJK s XJ*'YK 
FMN s XM*YN 

C EQUATION 5 

RMAsRHO»HS*B*PL/4,0 
C EQUATION 6 

HJK3ABS(1./<RMA*SQRT< (FNW< J,K)#FNW< J,K> -OMEGA* OMEGA )**2 
1+<2.*CI*FNW( J,K)*0MFGA>»*2) ) ) 

HMNsABS ( 1 ,/(RM ASSORT ( (FNW(M,N)#FNW(M,N)-0M£GA*0MEGA)**2 
l+<2,#CI*FNWfMiN)»0MFGA)**2) ) > 

C CALCULATE .JOINT ACCEPTANCES 

C EQUATION 7 

ACs ( 1 , ~ ( OMEGA/FNWI J » K ) >**2) Ml , - ( OMEGA/FNW ( M , N ) ) **2 ) *4 , *C I *C I *OMEG 
1A**2/(FNW( J,K)#FNW(M,Nn 
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8Cs-2.*( (C!*QMEGA/FNW< J,K> ) * ( 1 , -< OMEGA /FNW ( M , N ) ) *»2 > -C I *0 MEGa/ 
iFN W( M, OMEGA /FNW< J»K) )**?) ) 

CC=AC*AC+BC*BC 

CONI = (Al»OMFGA*PL/C)*( Al» OMEGA #Pt./C) 

CONSsR Ja-P I 
CON3*RM*PI 

CON6=COM2*CON3*( ( 2 , +EXP ( -Al*PL*OMEGA/C ) « { (-1. >**( J+l>*( <-l, )*# 
t<M + l> ) ) )/( <C0N1 + C0NP#C0N2)MC0N1 + C0N3*C0N3) ) ) 

CON7=RK*PI 

C0N8sRN*PI 

3 F (J ,NE. M) GO TO 46 

PJJM s {Al#PL* OMEGA/C/2, *(1,/(C ON 1+C0N2*C0N2)+ 
11,/(C0N1+C0N3*C0N3) )+C0N6) 

GO TO 47 

46 PJJM s C0M6+RM /(CON1+CON3«CON3) »< ( (-1, >**(J-M)-1. )/(?.* 
KRJ-RM) ) + ( (-1 , )**( J+M>-1 , ) / ( 2 . * { R J*RM 5 } ) 

2 +R J / ( CON1+CON2#CON2 )* ( C (-1 . ) **(M-J)-1 , ) / ( ? , a (RM-RJ) ) 

3*< (-1. )**<M*J>-1 . )/<2,*(RM+RJ) ) ) 

47 CONI =( A2*0MEGA*R /C ) * ( A 2#0MFG A*R / C ) 

C0N6=C0N7*C0N8« C(2.*6XPC-A2*8 *OMEGA/C>*<<-l.)**(K+l)+(<-l,)»* 
l(H+l>)))/( (CONl+CON7*CON7)*{CnNl+CON8aCONR) ) ) 

I F ( K ,NE. N) GO TO 48 

PJKN s ( A?*8 »0MEGA/C/2,a(1 ,/(C0Nl+C0N7aC0N7)* 

11 . /(C0N1+C0N8*C0N8 ) ) *C0N6 ) 

GO TO 23 

48 PJKN g C0N6+RN / < C0N1+C0N8*CQN8 ) * { ( ( -j , ) ** ( K-N ) -1 , ) / ( 2 . * 

1 (RK-RN) > + ( (-1 . . )/<2 .MRK+RN) > ) 

2 *RK /( C0N1+C0N7*C0N7)*< ( (-1 . )«#(N-K)-1 , >/<2 ,*(RN-RK)> 
3+<<-l,)*»(N + K)-l. ) /(2 ,* (RN + RK ) ) ) 

23 P0MEG=0MEGA/<2,*PI > 

P0WJ2=PJJM«PJKN*AC/9QRT{CC) 

PI WWsPJ WW + FJK*FMN*HJK#HMN*P0WJ2 
IF (POMEG .GT, 5, ) GO TO 22 
WRITE (3) POMEG , J , K , M , N , POW J2 

ICTslCT+1 
?? CONTINUE 

OMEG { I >*0MEGA/(?,*PI ) 

C calculate DISPLACEMENT, STRESS, and acceleration 

C EQUATIONS 8» 13, 14, 16, AND 17 

PIWW=PIWW*S»S*PIPP 
SWW( I )sPIWW*2.*PI 
PIWG =0MEGA**4*P I WW 
PIWGK I)=4.215093E-05*P!WG 
PSSWsGAM2*PlWW 
SSSF ( I )=2.*PI*PSSW 
C EQUATION 12 A 

SPPPt I ) =2 . *P I *P I PP 
IP ( IC ,NE. 1 ) GO TO 41 
C MEAN SQUARE EQUATIONS 

FlsOMEG(l) 

G1 = SWW< I > 

G2=SSSF( I) 

G3sP IWG1 ( I > 

G4sSPPP( I ) 

ATS1=0,0 
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ATS2*0 . 0 
ATS3sO , 0 
ATSAsO ,0 
IC = 2 

GO TO 21 

41 ATS1=ATS1*(G1+SWW< I > ) /2 . * < OMEG < I > -Fi ) 
ATS2sATS2+(G2+SSSF{ ! )1/2,*(0MEG( 

ATS3*ATS3+<G3+PIWG1< I ) ) /2 , * < OMEfi ( 1 )-Fl) 

ATS4 = ATS4+(G4 + SPPP t I ) ) /2 . * ( OMtG { I ) -FI ) 

C ROOT-MEAN SQUARE EQUATIONS 

FlsOMEG ( I ) 

GlsSWW { I ) 

G2sSSSF( !> 

G3sPlWGl( I ) 

G4 sSPPP ( I ) 

I F ( SWW ( I ) ,LT. SWMI ) SWMlsSWW(I) 

IF(SWW(I) ,GT. SWMX) SWMX*SWW<n 
IF(SSSFU) .LT, SSMI)SSMIbSSSF(I ) 

IF{SSSF( I) ,GT, S$MX)SSMX=SSSF( I ) 

I F { P I WG1 ( I > , LT . P l M I ) PIMIsPIWGlfn 
I F { P I WG1 ( I ) ,GT: PIMX) PIMX=PIWGl(n 

DSW a SWMX/1.E9 
DSS = SSMX/1.E9 
DPI a PIMX/1.E9 
OHEGAs OMEGA* CON 

21 CONTINUE 

KMs I 

GO TO 98 
99 KMs I -1 

C EQUATtON 12B 

98 PAS=ATS4 

ATlsSORT(ATSl) 

AT2sSQRT(ATS2) 

AT3sS0RT(ATS3> 

AT4sSQRT { ATS4) 

PA« A T4 

RUAsl70 .576+10 ,*A10G10( PAS) 

END FILE 3 
REWIND 3 
CALL PRNT 

I F ( DSW .LE. SWMI ) GO TO 302 
DO 301 1=1, KM 

IF(SWW<I> , LT i DSW) SWW(I)sDSW 

301 CONTINUE 

302 I F ( DSS .LE. SSMI) GO TO 304 
DO 303 I si # KM 

IF(SSSFCI) . LT , DSS) SSSF(I)sDSS 

303 CONTINUE 

304 I F ( DP I .LE, PIM!) GO TO 306 
DO 305 I si , KM 

I F { P I WG1 ( I ) ,LT. DPI) PTWGKDsDPI 

305 CONTINUE 

C PLOT EXCITATION, D I SPLACEMENT , STRESS, AND ACCELERATION 

306 CALL SCOUTV 

CALL 0UIK3V(-1,44,BCDX,BCDY1,-N3R,F3RD,SPPF) 
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oooooooooooooooooooo 


WRITE ( 16 » 109 ) RL A 
CALL SMXYV(l.l) 

CALL, QUlK3V(-l,44,BCDX*BCDY?»-KM,0Mf:G>SWW) 

WRITE (16,106) ATI, FINN 

CALL 0UIK3V<-l»44,BCDX,8CDY3,-KM,OMEG,SSSF> 

WRITE (16,106) A T2 , F J NN 

CALL GUI-K3V<-1,44,BCDX,BCDY,"KM,QMEG,P1WG1> 

WRITE ( 16,3 06 ) A T3 , F I NN 

CALL OU1K3V(-1,44,8CDX,BCDY4,-KM,OM6G,SPPP) 

WRITE (16,107) PA 
CALL CLFAN 

106 F0RMAT(11X,31H ROOT-MEAN-$OUARE RESPONSE = E15.7,5X,6HFTNM = 

1F6.?,10X,6HRFRPC1 ) 

107 FORMAT ( 11X , 30H ROOT-MEAN-SGUARE PRESSURE *,E15>.7) 

104 FORMAT(94HiDYNAMlC RESPONSE OF FOUR-SIDE FIXED RECTANGULAR SHELL P 
1ANELS CROSS-REINFORCED WITH STIFFENERS) 

109 F0RMAT(11X,24H0VERALL PRESSURE LEVEL * El5 , 7 , 8HQEC 1 PELS ) 

RETURN 

end 

SORInlN A 

SIBFTC RSF1 

SUBROUTINE RSF 
PROGRAM NUMBER B23-1002-7A 

pgp Pp | 

DYNAMIC RESPONSE OF TWO-OPPOSITE-SIDE SIMPLY-SUPPORTED 
AND OTHER TWO SIDES FIXED RECTANGULAR SHELL PANELS 
UNDER RANDOM PRESSURE FIELD 
INPUT PARAMETERS 
PL s PANEL LENGTH 
B * PANEL WIDTH 
RHO = MASS DENSITY OF PANEL 
HS = PANEL THICKNESS 
Cl s DAMPING RATIO 
X = COORDINATE OF VECTOR R 
Y = COORDINATE OF VECTOR R 

FINN = CONSTANT T-0 DETERMINE OCTAVE BAND CENTER FREQUENCIES 
A1 = CONSTANT 
C s SPEED OF SOUND 

BP = LENGTH OF PANEL UNDER EXCITATION 
PLP = WIDTH OF PANEL UNDER EXCITATION 
F3RD = ONE-THIRD OCTAVE BAND CENTER FREQUENCY 
S3RD = SOUND PRESSURE LEVEL 

DIMENSION FNW(10,10) ,F3RD(40) , S3RD ( 40 ) , BCDX ( 12 ) , 8CDY ( 1 2 ) 

DIMENSION F!PW(40> ,SPPF(40) ,OMEG(400) , P I WGl ( 400 ) , BC0Y1 (1?) 
DIMENSION 0CDY2(12) ,SWW(400) , SSSF(400) ,8 CDy 3(12) ,eCDY4 (12) 
DIMENSION SPPP (400) 

COMMON /OUTPUT/ FNW , F3RD , S3RD , F IP W , SPPF , OMEG , P l WGl , SWW , 
1SSSF,SPPP,N3R,KM, JK, I CT , PAS , ATS1 , ATS2 , ATS3 , A T1 , AT2 , AT3 , 
2QX,0Y.QW,GAM2, IRD 

COMMON /INPUT/ PL,8,RH0,HS,CI,X,Y,FINN,A1,C,PLP»BP, 
1A2,E,EP,VIP,AL1,BL1, All, AI2,H2,Hp,RAD,RH0P,RLA 
DATA 8CDX/72HFREQUENCY (HZ) 

1 / 

DATA BCDY/72HACCELERATION G SO/HZ 
1 / 
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DATA BCDY2/72HDISPLACFMENT INCH SQ/NERTZ 
i / 

DATA BCDY3/72MSTRESS spectral DENSITY psi sq/hfrtz 
1 / 

DATA BCim/72HSPPF DR/HZ 
1 / 

DATA BCDY4/72HSPECTRAL DENSITY psi SO/HERTZ 
1 _ / 

WRITE (6,108) 

SUM Is 1.E10 
SWMXs l.E-10 
SSM ! = 1 . El 0 
SSMXs l.E-10 
PIHIs 1.E10 
PIMXs l.E-10 
REWIND 3 
JK = 2 
ICT = 0 

PIs3. 1415927 
CALL CAMRAV<9) 

CALL SMXYV (1,0) 

C CALCULATE EXCITATION SPECTRAL DENSITY 

C EQUATIONS 11 AND 12 

DO 9 Iol,N3R 

SPPF< I )=S3RD( I )-10.*ALOG10(0,23l5*F3RD( I ) ) 

FIPW< I >*1, /<2.0*PI)*10 , 0**( ( SPPF { I ) -170 . 576 ) /1 0 ♦ 0 > 
CONTINUE 

CALCULATE NATURAL FREQUENCIES 
EQUATIONS 1 AND 2 
HCsE*HS*H$*HS/n2.#<l,-VlP*VlP) ) 

DYsHC+EP*AI?/ALl 
DXsHC+EP*AI 1/8L1 
BN s RHO*HS*RHOP*HP 
FCON s P I *P I /SORT ( BM ) 

DO 20 J=1 , 9 
DO 20 K = 1 » 9 
RJ = J 
RK s K 


I Ft J 

.NE, 

1 

.AND, 

K 

, NE . 

1) 

GO 

TO 

18 ' 

I F ( J 

.EG), 

1 

, AND, 

K 

,NF. 

1) 

GO 

TO 

16 

I F ( J 

.NE. 

1 

.AND, 

K 

. EO. 

1) 

GO 

TO 

17 


FHZ * DX*(1. /PL>**4 + DY*<i.5056/8)**4 + 2,*HC*<l./PL>*U,/PL>* 

1{1.1165/B)*(1.116'5/B}+E*HS/(RAD*rAD*{PI )*#4* (1 , + (1 . /1 . 5056>*(i , / 
21,5056)*(PL/S)#fPL/B) )**2) 

GO TO 19 

36 FHZ = DX«(1, /PL>**4+DYM (RK+,5)/B)##4 

1+2 . *HC*(1, /PL>**2 #( t (RK+,5)#< (RK+ , 5 ) - . 6366 ) ) /< B*B > 5 
2+E*HS/(RAD*RAD*(p! )*M»<l. + <l./(RK+.5}> **2»(PL/R)*(PL/3))*»2) 

GO TO 19 

17 FHZ = DX*( RJ /PL)**4+DY»(1,5056/B 5**4 
1*2, *HC* ( RJ/PL)*(RJ/PL )* (1 .1165/0) *( 1.1165/8) 

2+E*HS/(RAD*RAD*(P! )**4*<1,+(RJ/1,5056 > **2* < PL/B } * ( PL/B ) >**2> 

GO TO 19 

18 FHZ = DX»{ RJ /PL)**4+DY#( <RK+,5)/B)**4 
i+2,«HC*< RJ/PL)*(RJ/Pl)*( (RK*.5)*< < RK+ . 5) - . 6366) 
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2 +E*HS/(RA0#RAD»(Pn**4Ml. + < RJ /(RK+.5>)**2 
3*( PL/8)MPL/B) )**2> 

19 FNW { J , K ) = FCON*SQRT(FHZ> 

20 CONTINUE 

C CALCULATE CONSTANT TO CHANGE DISPLACEMENT INTO STRESS 

C EQUATION 16 

QX = 0 ,0 
QY = 0 .0 
QW = 0 ,0 
S=8P#PLP 
DO 49 IJ=1,5,2 
DO 49 IK=1,5,2 
RMsIJ 
R N = I K 

FSIN=SIN(RM*PI*X/PL>*SIN(RN*PI*Y/R) 

QCON' = RM*RN#(DX* (RM/PL )**4 + 2 ♦ «HC* ( RM»RN/ ( PL*B J) **2+DY« ( RN/B ) «*4 > 
QX=OX+( tRM*PI/pL>**? + VlP*<RN#PI/e>*»2>«(FSIN/OC0N) 

QY=GY+< <RN*PI/B>**2+VIPMRM*PI/PL>*»2)*(FSIN/QC0N) 

49 QWsOW+FS I N /QCON 
HI sHS+H2 

GAM2sF*H1*H1*E/(4,*(1.-VIP#VIP)**2)*( ( QV »QX+GY*Q Y ) / ( GW* GW ) ) 

C CALCULATE EXCITATION PRESSURE FROM ONE-THIRD OCTAVE DATA 

lOOfl CONs 2, 0/FINN) 

QMEGA=F3RD(1)*2.0*P! 

ICsl 

DO 21 1=1,400 
N3R1=N3R-1 
DO 40 JAsl,N3Rl 

IF ( OMEGA ,GT. ( F3RD ( N3R ) *2 , *P I )) GO TO 99 
IF(OMEGA-(F3RD(JA)*2,*P!))25»25,30 
25 KAsJA 

GO TO 35 
30 IAsJA+1 

IF(OMEGA-(F3RD(IA)*?,#PI )) 25, 25, 40 

35 PIPP = FIPW(KA> + CFIP^|(KA*1)-FIPW(KA) )»{0MEGA-F3RD(KA)*2,*PI )/< (F3RD 

1(KA+1)-F3RD<KA) )*2,*PI) 

GO TO 45 
40 CONTINUE 

45 PIWWsQ.O 
PIWWI=0.0 
C EQUATION 3 

DO 22 J = 1 » 5 
DO 22 K = 1 , 5 
DO 22 M = 1 , 5 
DO 22 N = 1 » 5 
RJ = J 
RKsK 
RMsM 
RN = N 

C CALCULATE THE NORMAL MODE 

C EQUATION 4 

XJsSlN(RU*PI*X/PL) 

XMbSIN(RM»PI*X/PL) 

CY1 s 1.5056#PI*Y/8 
CYK = (RK*.5)#PI*Y/R 
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CyN = (RN+,5)*PI*Y/B 

CHYJ s <EXP<CY1)+EXP(-CY1> >/2. 

SHY! s (EXP(CYU-exP(“CYl) )/2. 

CMYK = (EXp(CYK)+6XP(-CYK) )/2. 

SHYK s (EXP(CYK)-EXP(-CYK) >/2. 

CHYN = (EXP(CYN)+EXP(-CYN>)/2. 

SHYN = (EXP(CYN)~EXP(~CYN) )/2. 

211 IF (K .GT, 1) GO TO 212 

YK = CHYl^C0S(CYU-.9825*(SHYj -SIM(CYI) ) 

GO TO 215 

21? Y* = CHYK-COSCCYKl-(SHYK-SIN(CYK> > 

215 I F { N .GT. 1) GO TO ?16 

YN s CHY1-C0S(CY1>-.98?5*(SHY1-SIM(CY1) ) 

GO TO 217 

235 YN = CHYN-CnS{CYN)-(SHYN-SIN(CYN) ) 

217 FJK = XJ*YK 
FMN a XM*YN 
C EQUATION 5 

RMA»RH0*HS*B*PL/4 , 0 
C EQUATION 6 

HJKaABSf 1,/<RM ASSORT ( ( FNW ( J , K ) *FNW ( J , K ) -OMEGA IOMEGA ) **2 
1+(2.#CI*FNW< J,K)*0MEGA}**2) ) ) 

HMN=ABS< 1 * / t RMA*SORT ( (FNW{M*N)ttFNW(M,N)-C)MEGA*0MEGA)#«2 
1+<2,«CI»FNW(M,N)*0MEGA>**2) ) } 

C CALCULATE JOINT acceptances 

C EQUATION 7 

ACsf 1 ,-<OMEGA/FMW< J,K> . -<OMEGA/FNW(M, N)) **2>+4 , *C!»CI*QMEG 

1A»*2/(FNW( J,K)*FNW(M,N) ) 

BC=-2.*( (CI*OMEGA/FNW( J,K) >*< 1 OMEGA /FNW (M , N ) 3**2 ) -C I *OMEGA/ 

1FN W < M, N> *( 1.- COMEG A/FNWC J, K) )**2) ) 

CCs AC»AC+8C*BC 

CON3s(Al»OMFGA#PL/C)*( ai*omfga*pl/c ) 

CON2=RJ#Pl 

C0N3=RM»PI 

C0N6sC0N2*C0N3*( {?.*EXP(->A1*PL*0MEGA/C)*( ( *1 , ) **{ J + l ) + < < -1 , )** 
1CM+3 > ) ) >/C (C0N1+C0N2*C0N2)*(C0N1+C0N3*C0N3) ) ) 

CON7sRK«P! 

C0N8aRN*P I 

IF (J . NE . M) GO TO 46 

PJJM = (Al#PL»0MEaA/C/2.«(l,/(C0Nl*C0N2*C0N2)* 

li./CCONl+CON 3 *CON 3 ) J+C 0 N 6 ) 

GO TO 47 

46 PJJM s C0N6 + RM /( CONl + CON3*CON3 ) *C (( -1 .)** C J-M ) -1 .)/(?, * 
l(RJ-RM)) + ((-:l,)**(J«-M)~l.)/(2,#(RJ*RM))) 

2 +RJ /(C0N1+CQN2*C0N2)*( C (-1, )*»CM-J)-1, WC2,*CRM-RJ) ) 
3+((-1,)**(M+J>-1,)/C2.*(RM+RJ))) 

47 CONI =(A2*OMEGA*R /Cl *( A2#0MEGA*R / C ) 

C0N6sCGM7*C0N8#< (2,*EXP(-A2*B *OMEGA/C)»((-l,)**(K+l)+(C-l,)»* 

1(N + 1> > > )/{ (CONI + C0N7«C0N7)* (CONI + C0N8*C0N8)) ) 

IFCK , NE , N) GO TO 48 

PJKN = <A2*B *OMEGA/C/2,*(i f /{CONl*CON7*CON7>* 

11 ./(CON1 + O.OM8*CON8> )*00N6> 

GO TO 23 

48 PJKN = CGN6 + RN / ( C0N1 + C0N8*C0N8 1 * ( ( ( -1 . ) ** ( K-N ) -1 . ) / ( 2 , * 
l(RK-RN) )+( (-1, )**(K*N)-1, )/(2.*(RK*RN) > ) 
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2 *RK /<C0Nl+C0N7ftC0N7)*( (( -1 , )**(N-K 5-1 ,)/(?. #<RN-RK) > 

3 + M-l. ) ** { N±K ) -1 . ) 7 < 2 , * < RN + RK ) ) ) 

23 P0MEG=CMEGA/{2.*PI ) 

P0WJ2=PJJM*PJKN*AC/RQRT< CC) 
PIWWspiww*FJK*FMN*HJK*HMN*P0WJ2 
IF (POMFG . GT . 5. > GO TO 22 
WRITE < 3 ) POHEG , J , K , M , N , POW J2 

ICT^ICT+I 
22 CONTINUE 

OMEG ( I )=OMgGA/(?,*Pl ) 

C CALCULATE DISPLACEMENT, STRFSS , AND ACCELERATION 

C EQUATIONS 8 > 13 , 3 4 , 3 6, AND 17 

PIWWsPJWW*S*S*PIPP 
SWWU)sPIWW*2.*PI 
PIWG sQMEGA**4*PIWW 
PIWGK n=4.2l5 0R3F-05*plWG 
PSSNsGAM2*PlWW 

sssm )=2 ,*pi»pssw 

SPPP( I 5 =2 , *P I «P I PP 
IF ( IC .NE. 1 ) GO TO 41 
C MFAN SQUARE EQUATIONS 

F l = OM6G ( I ) 

GlsSWW < I ) 

G2 = SSSF( I) 

G3 = PIWG1( I ) 

G4 = SPPP ( I ) 

ATSlsO ,0 
ATS2s0.0 
ATS3sO , 0 
ATS4sO.O 
ICs? 

GO TO 21 

41 ATSlsATSl+(Gl + SWWU ) )/2.*(0M£G( I)-F1) 

ATS2aATS2+(G2 + SSSFM ))/?,*( OMEG M >-Fl) 

ATS3sATS3+(G3 + PTWG1 ( I) > 72 . * (OMEG < I) ~F1 ) 

ATS4sATS4+(G4 + SPPP ( I) ) 72 , * ( OMEG (I ) -FI ) 

C ROOT-MEAN SQUARE EQUATIONS 

F l = OMEG < I ) 

GlsSWW ( I ) 

GpsSSSFf I) 

G3sP[ WG1 ( ! ) 

G4sSPPP( I ) 

I F ( SWW ( I ) , LT , SWM I ) SWMI=SWW(I5 

IF(SWW(I) , GT , SWMX> SWMX=SWW<I> 

IF(SSSFM) .LT, SSMI)SSMI=SSSF(n 
I F ( SSSF < I ) .GT, SSMX>S$MX=SSSF( n 
IF(PIWGKI) .LT. P I M I ) P I M I =P I WG1 ( I ) 

I F ( P I WG1 ( I ) .GT, P I MX 5 p IMX«P I WG1 t I > 

OMEGA= OMEGA* CON 

21 CONTINUE 
KMsl 

GO TO 98 
99 KM=t-l 
98 PASsATS4 

ATlsSQRT(ATSl) 
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AT2=S0RT(ATS2> 

AT3=SQRT(ATS3> 

AT4=SGRT(ATS4) 

PAs AT4 

RLA=170 .576+10 ,*ALOG10( PAS) 

END FUF 3 
REWIND 3 
CALL PRNT 
DSW = SWMX/1,69 
DSS s SSMX/1.89 
DPI s PIMX/1.E9 
DO 301 1=1, KM 

IF(SWW(I) , LT , DSW) SWW(1)=DSW 
I F < SSSF ( I ) .LT, DSS) SSSF(I)sDSS 
IF(PIWGKI) ,LT . DPI) P!WGlM)snpI 
301 CONTINUE 

C PLOT EXCITATION, DISPLACEMENT, STRESS, AND ACCELERATION 

CALL SCOUTV 

CALL QUIK3VC'*1»44i80DX,RCDY1>-N3R,F3RD,SPPF) 

WRITE M6.109JRLA 
CALL SMXYV(l.i) 

CALL 0UIK3V(-1,44,BCDX,RCDY?,-KM,OMEG,SWW) 

WRITE (16,106) ATI, FINN 

CALL QUIK3V(-1,4 4,BCDX,8CDY3,’-KM,0MEG,SSSF) 

WRITE (16,1065 AT2 , F I NN 

CALL 0UIK3V(-1,44,8CDX,8CDY,-KM,0MEG,PIWG1) 

WRITE (16,106) AT3 , F I NN 

CALL QUIK3V(-1,44,BCDX,BCDY4,-KM,0MEG,$PPP) 

WRITE (16,107) PA 
CALL CLEAN 

106 FORMAT (HX,3lH ROOT -MEAN-SQUARE RESPONSE = E15 , 7 , 5 X , 6HF I NN = 

1F6,2,10X,6HRSFRP1 ) 

107 F ORMAT ( 11X, 30H ROOT-MEAN-SQUARE PRESSURE =,El5,7) 

108 FORMAT(107H1DYNAMIC RESPONSE OF TWO-OPPOSITE-SIDE SIMPLY-SUPPORTED 
1 AND OTHER TWO SIDES FIXED RECTANGULAR SHFLL PANELS,/, 

228H UNDER RANDOM PRESSURE FIELD) 

109 FORMAT (11X,24H0VERALL PRESSURE LEVEL = E15 , 7 , 8HDEC I 8ELS > 

RETURN 

END 

SDATA 
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I. METHODS OF VERIFICATION 


Comparison of the computed results with experimental data is the best 
verification. Verification can also be obtained by hand calculation of the 
responses according to the simple formulas by assuming the structure vibrates 
m its fundamental mode. 

The root-mean-square responses are computed, printed, and plotted out 
as output. Good engineering judgement on these rms responses may also serve 
as a check of the results. 
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SECTION III. DECK SET-UP 


A. 


Computer Configuration 

1. Computer 

2. Core Size 

3 . Language 

4. Operating System 

5. Plotter Required 

6 . Punch 

7 . Tape Assignments 


IBM 7094 
32K 

FORTRAN IV 
IBSYS 
SC4020 
NO 


Physical 

Unit 


Logical 

Unit 


System 

Function 


A2 5 

A4 3 

A8 

B1 6 

B3 2 


Input 

System Scratch 
SC 4020 Output 
Print Output 
System Scratch Ov erla y 


B. Estimated Running Time 

Execution time is approximately 55 minutes for all three boundary 

conditions, for frequency range 5 to 5000 Hertz, and frequency increment 
FINN =33. 

C. Restart Procedure 
None 

D. Deck Sequence 


1 . 

$JOB CARD 


2. 

$ EXE CUTE CARD 

3. 

$IBJOB CARD 


4. 

$IBFTC CARD 


5. 

SOURCE DECK 

(RANDOM) 

6. 

$IBFTC CARD 


7. 

SOURCE DECK 

(SUBROUTINE PRNT) 

8. 

$IBFTC CARD 


9. 

SOURCE DECK 

(SUBROUTINE GRIDIV) 
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10. $ORIGIN CARD 

11. $IBFTC CARD 

12 SOURCE DECK (SUBROUTINE RSR) 

13 $0RIGIN CARD 

14. $IBFTC CARD 

15. SOURCE DECK (SUBROUTINE RFR) 

16 . $ORIGIN CARD 

17 . $IBFTC CARD 

18. SOURCE DECK (SUBROUTINE RSF) 
19 $DATA CARD 

20. DATA DECK 
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E. INPUT DATA 


Refer to Figure 1 for geometric dimensions. 


NOTATION 

DEFINITION 


CARD 

FORMULAS MNEMONICS 


NO. 

FORMAT 



NP 

Program desired - NP = 1 (RSR) 
NP - 2 (RFR), NP = 3 (RSF) > 

NP = 4 (ALL) 

i 

15 


N3R 

Number of data points m the one- 
third octave band excitation 
spectrum 

i 

15 


PL 

Axial length of panel (inches) 
(Along x-axis) 

2 

E15.8 

b 

B 

Width of Panel (inches) 
(Along y-axis) 

2 

E15.8 

P 

RHO 

Mass density of panel skin 
(lbf-sec 2 /m ^) 

2 

E15.8 

h 

HS 

Thickness of panel skin (inches) 

2 

E15.8 


Cl 

Damping ratio of panel 

3 

E15.8 

X 

X 

T>- 

Coordmate of r (inches) 

3 

E15.8 

y 

Y 

Coordinate of ^ (inches) 

3 

E15.8 

n 

FINN 

One-nth octave frequency 
increment 

3 

E15.8 

A i 

Al 

Correlation decay constant in 
axial length-direction 

4 

E15.8 

c 

C 

Speed of sound (m/sec) 

4 

E15.8 

V 

PLP 

Length of panel subjected to 
excitation (inches ) 

4 

E15.8 

b' 

BP 

Width of panel subjected to 
excitation (inches) 

4 

E15.8 


COL. 

Col. 5 

4 

Col. 6-10 | 

Right Justi-* 
fied 

Col. 1-15 

Col. 16-30 

Col. 31-45 

Col 46-60 

Col 1-15 

Col. 16-30 
Col. 31-45 
Col. 46-60 

Col. 1-15 

i 

Col 16-30 
Col. 31-45 
Col. 46-60 
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E. INPUT DATA (Continued) 


NOTATION 

DEFINITION 



CARD 


FORMULAS 

MNEMONICS 


NO. 

FORMAT 

COL. 

A 2 

A2 

Correlation decay constant in 
circumferential width-direction 


5 

E15.8 

Col. 1-15 

E 

E 

Young’s modulus of panel skin 
(lbf/m^) 


5 

E15.8 

Col. 16-30 

E' 

EP 

Young T s modulus of stiffeners 
(lbf/m 2 ) 


5 

E15.8 

Col. 31-45 

V 

VIP 

Poisson’s ratio of panel skin 


5 

E15.8 

Col. 46-60 

a l 

ALl 

Spacing of width-direction 
stiffeners (inches) 


6 

E15.8 

Col. 1-15 

b l 

BLl 

Spacing of length-direction 
stiffeners (inches) 


6 

E15.8 

Col. 16-30 

*1 

All 

Moment of inertia of one length- 
direction stiffener with respect 
to neutral axis (xnch^) 


6 

E15.8 

Col. 31-45 

I 2 

AI2 

Moment of inertia of one width- 
direction stiffener with respect 
to neutral axis (inch^) 


6 

E15.8 

Col. 46-60 

h 2 

H2 

Largest height of stiffeners at 
point investigated (inches) 


7 

E15.8 

Col. 1-15 

h' 

HP 

Smeared- out thickness of stiffeners 
(inches) 

7 

E15.8 

Col. 16-30 

a 

RAD 

Radius of panel (inches) 


7 

E15.8 

Col. 31-45 

P 

RHOP 

Mass density of stiffeners 
(lbf-sec 2 /in^) 


7 

E13.8 

Col. 46-60 

f 

F3RD 

Frequencies (Hertz) 

4 Values 
Per Card 

F8.1 

Col. 1-8,17-24 
33-40, 49-56 - 

s 3 r (i) 

S3RD 

One- third octave pressure level 4 Values 
spectrum of excitation (decibels)Per Card 

F8.1 

Col. 9-16,25-32 
41-48, 57-64 [ 



Use as many cards as necessary for 
F3RD and S3RD 
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F. Restrictions and Limitations 

N3R <_ 40 
FINN <_ 38 

G. Diagnostics 
None 

H. Quanity of Output 

For case with frequency range of 5-5000 Hertz and frequency increment 
FINN = 33, the printed output will be 20 pages per boundary condition. Plots will 
be 5 per boundary condition. 
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I* 


Output Definitions 


NOTATION 


FORMULAS 

ft 

b 

P 

h 

^jk 

x 

y 

n 

A i 

c 

Z’ 

b' 

A 2 

E 

E' 

v 

a l 

H 

*1 


MNEMONICS 

PL 

B 

RHO 

HS 

Cl 

X 

Y 

FINN 

Al 

C 

PLP 

BP 

A2 

E 

EP 

VIP 

AL1 

BL1 

All 


DEFINITION 


Axial length of panel (inches) 

Width of panel (inches) 

Mass density of panel skin 

(lbf-secVm 4) 

Thickness of panel skin (inches) 

Damping ratio of panel 
Coordinate of r (inches) 

4 - 

Coordmate of r (inches) 

One-nth octave frequency increment 

Correlation decay constant m axial length- 
direction 

Speed of sound (in/ sec) 

Length of panel subjected to excitation 
Cinches) 

Width of panel subjected to excitation 
(inches) 

Correlation decay constant m circumferential 
width-direction 

2 

Young’s modulus of panel skin (lbf/m ) 

Young’s modulus of stiffeners 
(lbf/in ) 

Poisson's ratio of panel skin 

Spacing of width-direction stiffeners (inches) 

Spacing of length-direction stiffeners (inches) 

Moment of inertia of one length-direction with 
respect to neutral axis (inch^) 
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I. Output Definitions (Continued) 


NOTATION 


DEFINITION 

FORMULAS MNEMONICS 


*2 

A12 

Moment of inertia of one width-direction 
stiffener with respect to neutral axis 

(mch^) 

h 2 

H2 

Largest height of stiffeners at point 
investigated (inches) 

h' 

HP 

Smeared-out thickness of stiffeners (inches) 

a 

RAD 

Radius of panel (inches) 

p 1 

RHOP 

Mass density of stiffeners 
(Ibf-sec^/m. 4) 

f 

F3RD 

Frequencies (Hertz) 

s 3l m 

S3RD 

One-third octave pressure level spectrum 
of excitation (decibels) 

s pp ( f > 

SPPF(I) 

Excitation spectral density (db/Hz) 

♦ pp <“) 

FIPW(I) 

Excitation spectral density (psi /rad/sec) 

j 

J 

Mode 

k 

K 

Point 

f 

FHZ 

Natural frequencies (Hz) 

w jk 

FNW ( J , K) 

Natural frequencies (rad/sec) 

f 

OMEG(I) 

Frequency (independent variable) (Hz) 


SWW(I) 

Displacement spectral density (inch /Hz) 

s ao^» f ) 

SSSF(I) 

Stress spectral density (psi /Hz) 

®ww^ r ’^ 

PIWGl(I) 

Acceleration spectral density (G 2 /Hz) 

s pp (f) 

SPPP(I) 

Excitation spectral density (psi z /rad/sec) 

W 2 (?) 

ATS1 

Mean square displacement (men ) 

o 2 (r) 

ATS 2 

2 

Mean square stress (psi) 

2 "*■ 
» (r) 

ATS 3 

2 

Mean square acceleration (g ) 

-»■ 

w(r) 

ATI 

Root-mean square displacment (inch) 
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I 


I. Output Definitions (Continued) 


NOTATION 


DEFINITION 

FORMULAS 

MNEMONICS 


o(r) 

AT2 

Root-mean square stress (psi) 

G(r) 

AT3 

Root-mean square acceleration (g) 

Q x 

QX 

Quantity for the calculation of y 2 

Qy 

QY 

2 

Quantity for the calculation of y 


QW 

Quantity for the calculation of y^ 

9 

y 2 (r) 

GAM2 

Constant to change displacement spectral 
density into stress 

J ,k,m,n 

J,K,M,N 

Mode Indices 

t 2 

J jkmn 

POWJ2 

Joint acceptance squared 
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J. Operator Ins true txon Card 


7094 - 


INSTRUCTIONS 


NAME 


Moore 


OP CODE 


Mu 


STACK T 


BIN 


^2lfe [go 4^ A? 


JOB 


IF EXCEEDS MAX 


ms3o 


FAST TAPES A B C D 


□ STR OST£ OOMPQRETSY 

INPUT TAPES 

1 

WORK 

LOGIC 


1 LOGIC I 

REEL NO 

ESI 


JfgOMPL /ASSMBL 
Jl£XECUTE 
Q PUNCH IBCD BIN) 









V4FTRN 

xJeftrn 

a aft 

□ PERT 

OMAP 

OFAP 

□ SCAT 

□ OTHER 










1 


i 



LINES OF OUTPUT OOOO's) 

DQ-i Q 6 - 13 O 15-30 DQVERl 


MAXIMUM TIME 

Lours minutes 


SS. 


programmer comments 


NUMBER OF CASES „ 


OVER 



□ SEE ON-LINE 


□ SEE TECHNIQUES 


□ MAX EXCEEDED 

OPERATOR COMMENTS 

D RETURN TO SYS 


D LINE MAX 


OPER IMIT 
OVER 


OUTPUT TAPES ONLY 

4020 

REEL NO 

! LOGIC 

BO 

UNIT 

NO OF CPYS 

SAVE 

TAPE 


B-l 

11 






rwm 





V 






I 


















■ 

■ 





NO FILES 

NO FRAMES 




COPY FLO | 

KALVAR 

1 

J5„ 

a 

a 

□ 

o 

o 

W'M 

F 




! .. 1 

i i 

i i 

I j 


1 



MSFC “ Form 533 fRev February 1966) 
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K. Save Labels 
None 
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INPUT FORM (SAMPLE CASE) 


10 


13 


20 


25 


30 


33 


40 


43 


30 


33 


SO 


03 


70 72| 


SAMPLE OUTPUT 



SAMPLE OUTPUT OF PROGRAM RANDOM 


FREQUENCY DISPLACEMENT STRESS- ACCELERATION ExCiTiliON 

RESPONSE RESPONSE RESPONSE 

HERTZ INCH SQ/HERTZ pS-I SQ/HERTz 6 SQ/HERTZ PSI SQ/RAD/SEc. 

5 *00000 0 • 2Z 1 06 1 7 IE-06 n* 33 **02308E 03 0* J 4445757E-05 0 *6(j 1 006Z6E-04 

5*00000 0*22S0617|E-06 0.33402308E 03 0, 1 4445757E-Q5 0 • 6q l 0U626E-Q4 

5*1 1687 0.22H97937E-06 0.33994264E 03 0. I 6 1 25262E-05 Q.6 J 4938U2E-04 

5 • 2 3 647 Q.22-894313E-06 0.3 4-5-9 3-1-8 6E 03 0. 1 7998 1 98E-Q5 0*629 1 9542E-04 

5.35887 0.23295265E-06 0.35199022E 03 0. 20086596E-05 Q.6<i 3 7a6U8E-Q4 

5 *4 a 4 1 2 0*237007*»6E-06 0.3581 1703E 03 0. 224 1 4964E-Q5 0 *6587 1 7 7 8E-Q4 

0^4U0694E-06 0.36431 132E 03 0 * 250 1 Q540E-Q5 Q.67399849E-Q4 

50*39679 0 • 2 1 680200E-06 0.32758667E 03 0* 1 4622498E-Q 1 0* 1 90D4992E-03 

51 .57476 0.-2-1 3596H32E’-06 0.32274247E CT3 0 • 1 580 1 1 52E-0 1 0 . 1 9003523E-03 

52*78026 0.21047653E-Q6 0.31802892E 03 G * 1 707798 3E-Q 1 0 » 1 9U02020E-03 

54*01394 D.20744474E-06 0.31344790E 03 0 • l 8 4 6 1 7 4 1 E -0 1 0 • 1 900Q482E-03 

55*27646 0 • 20450 1 90E-Q6 0.309Q0129E 03 0* 1 9962038E-0 1 0. 1 8 998907E-Q3 

56 *56848 0 . 20 1 64937E-06 0.3Q469H4E 03 0. 2 1 589457E-0 J 0* 1 8997296E-Q3 

57*89071 0* 19888860E-0A 0.30051962E 03 0 * 23355655E-0 1 0 * 1 8995647E**03 

59 *2*1384 0 • 1 9622 1 06E-06 0.29648898E 03 0* 25273482E-0 I 0 • 1 8 9 93960E-03 

60*62859 0 • 1 9364840E-06 0*29260l7l£ 03 0* 27357 1 4 1 E-0 1 0 • 1 8 9 92233E“03 

62*0*1572 0*191 1 7232E-06 0.28886037E 03 0 * 29622325E-0 1 0 * 1 8 990466E-03 

192.48369 0 . 642 1 2375E-06 0.970245ME 03 0.92159750E 01 0 ♦ 9 A834247E-04 

196*98278 0 . 88668852E-06 0.13397817E 04 0. 13958249E 02 Q • 9g6 l 7973E-04 

201*58703 0. 13158538E-05 0.19882482E 04 0.22719816E 02 0 • 94200486E-04 

206.29890 Q • 209 9 8 0 1 OE-05 G.31 727882E 04 0.39766065E 02 0 * 924 1 3703£-0<* 

21 1 *12091 0.33107326E-05 0.50024996E 04 0 f 68769506E 02 0 * 9 q585 1 55E-04 

216*05563 0 • 365955 1 6E-05 0.55295632E 04 0.8337S221E 02 0* 887 1 3867E-Q4 

221*10569 0* 2306 1 066E-QS 0.34845150E 04 0.57626974E 02 0 • 8 67 9 8839E-04 

226*27378 0 . 1 1 954749E-05 0.18063562E 04 0»32766058E 02 0 • 8^83905 1 E-04 

231 *56268 0 * 6465400 1 E-06 0.97691854E 03 0»J9436429E 02 0 • 82833453E**04 

236*97520 0 • 37a 34 1 90E-06 0.57167262E 03 0.12475065E 02 Q . 8 q780978E-04 

570*17398 0.42722072E-09 0.645S2825E 00 0*47209420E 00 0* 2Q496707E-04 

583*50116 0.49555338E-Q9 0.74877855E oq 0*600626aQE 00 0 • 1 9290879E-Q4 

597 • 1 398S 0.50793323E-09 0.76748443E oo 0»67523918E 00 Q* 1 8056867E-04’ 

611*09732 0* 32488827E-.Q9 0.49090447E 00 0.47372068E 00 0 * U7 940 1 0E*04 

625.-38104 0* 153-7-2795E-09 0 . 232282-I2E -00 0»2458546eE 00 0 * 1 S50 1 636E-04 

639*99863 0 . 946283 1 2E- 1 0 0.14298288E 00 0tl6599Q8lE 00 0 • 1 469 1 1 50E-Q4 

654*95788 0 « 844 1 8406E- 1 0 0.1275S577E 00 0»1624*l92lE 00 0* 1 #-1 Q3838E-04 

670*2*678 0.94 J98362E-10 0.14233323E 00 0«19878378E QO Q . 1 3SQ27B8E-04 

685.93353 0* 1 1 849360E-09 0.17904321E 00 0.27426463E 00 0 . 1 288769 1 E-04 

701*96645 0. 16238162E-Q9 0.24535777E 00 Q.41223900E QO 0 « 1 225821 6E-Q4 

4063*73444 0 . 33002897E- 1 6 0 . 4986720?E-0-7~ 0 » 94 1 02274E-04 0 • ft 57 7557E-06 

4158*71967 0 • 27798077E-* l 6 0 • 420Q2749E-07 0 t 86936097E-04 Q . 1 1 08 7 1 27E-06 

4255.925-1-1 0 . 2336-3096E--1 6 0 . 3530 IS-l 5E-07 0 * 801 40694E-Q4 -Q . 1 0585232E-06 

4355*40265 0 . 1 9584 1 42E- 1 6 0 , 29591 536E-07 0 » 736 82540E-04 0 . 1 Q07 1 606E-06 

4457*20532 0 • 1 6364966E- 1 6 0.24727378E-07 0 * 67532444E-Q4 0 • 9&459756E-07 

4561*38745 0 . 1 3623937E- 1 6 0 . 20585697E-07 0 . 6 1 66 479 2E-04 0 • 900805 B8E-07 

4668*00482 0 * 1 129-1 6 8~2E— 1-6 0 * 1 706 1 6-7 lE-07- 0 .-5605 7096E-04 0 * 8«(575686E-07 

4777* 1 1420 0.93091068E-17 0 . 1 40660 U E-07 0*506894 1 8E.-04 0* 7a942 1 1 2E-»07 

4888 *-7^39-3 0*^2Sft00-3E*»-l 7 -Q-. 1 1-5 2-2 &4 4 £-• Q 7- Qr4S5440-91-Ef04 0.-73 1 76862E"Q7 

MS VAL-UE 0. 1-2538-1-7-2E-1Q3 0 , 1 8^45 1 1-6E 06 Q.24598280E 04 0* 5 t 435825E-0 1 

RMS VALUE 0*11 197398E-01 0.43525988E 03 0.49596653E 02 Q.22679467E 00 

QX ■ 0.37250E-Q1 QY * 0.29523E-01 QW » 0.739Q3E 01 GAMMA2 « 0*l5ll0E 10 
FINN -«-30.Q 
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N . Sample Plots 


Figure 2 through 4 are three sample plots of this program. 
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FIGURE 1. GEOMETRY OF RECTANGULAR CYLINDRICAL SHELL PANEL 
CROSS-REINFORCED WITH STIFFENERS 






RUN 9582 

ROOT- ME AH- SQUARE RESPONSE “ G .11 19740E-01 


riHH s 30 00 


RffiPCl 


0 

1 

S 

P 

L 

A 

C 

t 

M 

C 

N 

T 


1 
N 
C 
H 

3 

a 

/ 

H 

E 

R 

T 

2 


i.OOXlo' 0 * 


1.00X10~ 0? 


1 .00X10 


1 .00X10 


-to 


1 .00X10 


l.OOxlO ** 11 


l.OOXto ” 11 


1 .00X10 


-IS 


1.00X10 


-14 



10.00 


100.00 

frequency chz> 


1000.00 


FIGURE 3. SAMPLE PLOT: DISPLACEMENT SPECTRAL DENSITY AT CENTER 

OF FOUR EDGES CLAMPED RECTANGULAR CURVED PANEL CROSS- 
REINFORCED WITH STIFFENERS 
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RUN 9582 

ROOT-KEAH-saUARE RESPONSE = 0 3QS2289E 03 FINN i 30.00 RSFRP1 



FIGURE 4. SAMPLE PLOT: STRESS SPECTRAL DENSITY AT CENTER OF 

TWO OPPOSITE EDGES SIMPLY-SUPPORTED WHILE OTHER TWO 
CLAMPED RECTANGULAR CURVED PANEL CROSS-REINFORCED 
WITH STIFFENERS 
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DISTRIBUTION 


This report is to be distributed as follows 


PR-SC 

MS-IL 

MS-T 

MS-I 

R-P&VE-RI 


1 

1 

1 

1 

12, plus reproducible master 
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